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ABSTRACT 

This  thesis  conducts  an  experimental  investigation  of 
the  effects  on  separation  of  the  closing  angle  of  the  stern 
section,  the  depth  of  submergence  or  proximity  of  the  free 
surface,  and  the  speed  of  a  submerged  body  of  revolution. 
The  results  of  this  investigation  are  then  compared  to  the 
results  obtained  from  two  existing  theoretical  criteria  for 
the  prediction  of  the  separation  point  on  such  bodies  in 
turbulent  flow.   The  point  of  separation  is  observed  experi- 
mentally by  photographing  dye-flow  patterns  over  the  stern 
of  the  model  at  three  speeds  from  .75  to  1.7  knots,  at  three 
depths  of  submergence  corresponding  to  surfaced,  near  sur- 
face and  fully  submerged  conditions  for  five  different  body 
stern  shapes  increasing  in  streamlining  from  a  blunt  hemis- 
phere to  a  finely  shaped  almost  pointed  section.   Equivalent 
body  profiles  to  the  outside  of  the  turbulent  boundary  layer 
displacement  thickness  are  computed  and  potential  flow  pres- 
sure-velocity distributions  based  upon  the  equivalent  pro- 
files are  calculated  and  used  in  evaluating  the  separation 
criteria. 

It  is  concluded  from  the  results  that  initiation  of 
separation  occurs  at  a  nearly  constant  body  angle  on  all 
bodies  tested  regardless  of  the  degree  of  streamlining.  Ex- 
perimental results  show  no  apparent  effect  of  depth  of  sub- 
mergence on  the  separation  point  and  the  effect  of  speed 
shows  no  clear  cut  pattern.   The  sensitivity  of  the  variables 
used  in  the  separation  criterion  on  the  predictions  of  the 
separation  point  are  evaluated  and  recommendations  are  made 
for  further  studies  and  experimentation. 
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NOTATION 

Symbol  Meaning  Units 

C  NOn- dimensional  coefficients 

E  Error  functions 

K  Matrix  functions  of  body 

K  Matrix  functions  of  prolate  spheroid 

L  Body  length  ft.  or  in. 

P  One  half  the  perimeter  of  the  body  profile   ft. or  in. 

R  Gaussian  weighting  functions 

Rs  Reynold's  number,  along  body  profile 

Rx  Reynold's  number,  along  x-axis 

S         Surface  area,  over  body  ft.^or  in. 2 

s 
U         Total  velocity  along  body  ft. /sec. 

Uoo  Free  Stream  velocity  ft. /sec, 

V  Velocity  of  body  ft. /sec. or  knots 
X  Distance  along  x-axis  ft. or  in. 

Y  Distance  along  y  -  axis  ft. or  in. 

ds  Element  of  arc  length  along  body 

g(x,t)     Function  of  prolate  spheroid  whose  ends  coincide 
with  body  ends  and  intersects  body  at  x  -  t 

g         Acceleration  due  to  gravity 

k         Matrix  functions  of  body 

k^         Longitudinal  virtual  mass  coefficient  of  given  body 

kf  (x)      Longitudinal  virtual  mass  coefficient  of  prolate 
spheroid  at  x  -  t 

n         Reciprocal  of  exponent  in  boundary  layer  velocity 
profile 

k'         Matrix  functions  of  prolate  spheroid 

viii 


NOTATION  (Cont'd) 
Symbol  Meaning  Units 

p         Dynamic  pressure  lb/ft2 

2 

q         Stagnation  pressure  lb/ft 

r         Radius  of  body  profile  (calculation     in. 
for  corrected  radius) 

r         Distance  from  location  of  point         in. 
source  on  x-axis  to  point  on  body 
profile  (calculation  of  pressure 
distribution) 

t         Location  of  point  source  on  x-axis,  also 

equal  to  position  of  Gaussian  abscissae  on 
x-axis,  normalized 

u  Component  of  velocity  along  x-axis      ft/sec 

v  Component  of  velocity  along  y-axis      ft/sec 

w  Axisymmetric  potential  function 

x  Distance  along  x-axis,  normalized 

y  Distance  along  y-axis,  normalized 

0  Angle  between  outward  normal  to  body 

and  x-axis  (calculation  of  pressure 
distribution) 

P  Constant  in  separation  criterion 

#  Angle  between  normal  to  body  and  y-axis  or 

between  tangent  to  body  and  x-axis 

A  Ratio  of  £*/s 

S  Boundary  layer  thickness  in. 

c* 

o  Boundary  layer  displacement  thickness      in. 

£  Boundary  layer  momentum  thickness  in. 

X         Ratio  of  body  length  to  body  diameter 

V  Kinematic  viscosity  of  water  ft  /sec 
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NOTATION  (Cont'd) 

Symbol  Meaning  Units 

£         Normalized  position  of  Gaussian  abscissae 

based  on  body  length  of  2  and  from  -1  to  +  1 

P  Mass  density  of  water  slugs/ft 

0         Velocity  potential  function 

T  Stokes  stream  function 

J I  Ratio  of  /£  in  calculation  of  corrected 

body  profile 

-O.         Stream  function  of  a  unit  source  on  the  x-axis 

Subscript  Meaning 

d         Subscripts  drag  coefficient 

f         Subscripts  coefficient  of  friction 

i         Gaussian  abscissae  index 

j  Gaussian  abscissae  index 

k         Velocity  is  in  knots. 

n         Refers  to  direction  of  body  normal  in  derivation 
of  pressure-velocity  distribution 

n  Index  of  iteration  for  solving  equation  of 

pressure- velocity  distribution 

p  Subscripts  pressure  coefficient 

pi  Quantity  referred  to  flat  plate 

r  Subscripts  residual  resistance  coefficient 

t  Subscripts  total  drag  coefficient 

x  Quantity  referred  to  x-axis 


I   -    INTRODUCTION 

Separation  of  turbulent  flow  with  its  associated  causes 
and  effects  has  for  the  last  half  century  been  a  subject  of 
controversy  in  the  field  of  naval  architecture  and  ship  de- 
sign.  It  is  in  this  aspect  of  fluid  flow  that  the  naval 
architect  has  his  greatest  lack  of  knowledge.   There  is  even 
considerable  debate  and  controversy  as  to  just  what  is  known 
and  what  is  not  known  on  the  subject  of  separation,  its 
causes  and  effects.   Little  has  been  mentioned  in  the  past 
regarding  this  aspect  of  fluid  flow  in  references  or  texts 
on  ship  design.   What  has  been  stated  has  been  conjecture 
based  on  limited  theoretical  calculations,  sketchy  experi- 
mentation, or  recalled  experience  and  observations. 

The  causes  and  effects  of  separation  and  the  para- 
meters relating  to  them  are  complicated  and  closely  linked 
and  must  be  fully  understood  before  any  definite  design  cri- 
teria of  value  can  be  established  for  application  to  ship 
and  submarine  forms.   Separation  of  flow  on  such  forms  in- 
volves the  three  dimensional  boundary  layer  and  turbulence 
theories  which  have  hardly  been  more  than  elaborated  be- 
cause of  their  associated  mathematical  difficulties.   Even 
the  mathematically  simpler  axialsymmetric   flow  becomes  in- 
volved with  the  problems  of  turbulent  flow  which  is  not  fully 
understood  or  theoretically  formulated  for  the  relatively 


simple  case  of  two  dimensional  flow  over  a  flat  plate. 

The  ultimate  goal  is  to  establish  a  valid,  but  simply 
applied,  theoretical  criterion  to  predict  the  separation 
point  and  a  foolproof  method  of  testing  such  a  criterion 
with  models  to  aid  naval  architects  in  improving  ship  de- 
sign.  To  do  the  latter  the  scaling  effects  involved  in 
separation  must  be  firmly  established  and  understood  as  well 
as  understanding  the  basic  phenomenon  of  separation  itself. 
The  most  important  and  difficult  step  in  realizing  the  for- 
mulation of  such  a  criterion  is  the  determination  of  the  sig- 

ta 
ni  fie  ant  parameters  to  be  used  and  the  quant  iAtive  effect  of 

the  location  of  the  point  of  separation. Faired  surfaces 
with  modest  closing  angfes  apparent ly  retard  separation  but  the 
limiting  magnitude  of  the  body  closing  angle  or  extent  of 
fairing  required  to  minimize  separation  is  unknown. 

The  effects  of  hydrostatic  pressure  and  atmospheric 
pressure  on  separation  in  general,  and  in  particular  the 
effect  on  the  limiting  closing  angle  have  been  the  subject  of 
considerable  debate.   Saunders  (Y)  gives  some  information  re- 
lating to  these  topics  but  this  information  has  apparently 
never  been  checked  or  verified.   The  effects  of  the  free  sur- 
face on  separation  on  a  submerged  body  as  it  approaches  and 
penetrates  the  free  surface  is  another  area  of  considerable 
debate.   Body  speed  and  surface  roughness   are  other  possible 
parameters. 


*  Numbers  in  brackets  refer  to  works  referenced  in  Appendix  H. 
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The  object  of  this  thesis  is  to  investigate  the  ef- 
fects on  separation  of  varying  depth  of  submergence  or  proximi- 
ty to  the  free  surface,  body  speed,  and  the  closing  angle  of 

i 
the  stern  of  the  body.   The  results  are  necessarily  qualitative 

in  nature  since  the  primary  tool  used  in  this  investigation, 
dye-stream  photography , is  limited  in  the  exactness  to  which 
the  point  of  separation  can  be  determined  by  noting  the  lo- 
cation of  the  familiar  flow  pattern  of  separation,  i.e. , 
abrupt  departure  of  the  streamlines  of  flow  from  the  body 
surface.   Originally  it  was  hoped  to  relate  the  observed 
movement  of  the  separation  point,  if  any,  with  the  independent 
variation  of  the  parameters  listed  above  to  the  form  drag  on 
the  body.   However,  the  results  of  the  drag  measurements  on 
the  model  were  unsatisfactory  due  to  faulty  or  inadequate 
instrumentation  and  this  could  not  be  done.   Two  separation 
criteria  ,  one  developed  by  Stratford  (2)    and  one  described 
by  Granville  in  (3)  ,  are  theoretically  evaluated  and  the  re- 
sults compared  to  the  observed  point  of  separation  from  the 
dye-stream  photographs.   This  provides  some  insight  as  to 
validity  of  the  criteria   and  the  significance  of  the  para- 
meters used  therein. 

It  is  felt  that  the  results  of  this  thesis  are  a  step, 
although  incremental,  toward  determining  what  the  significant 
parameters  in  separation  of  flow  are  and  how  they  are  related 
to  the  phenomenon  of  separation. 


II    -   PROCEDURE 

A.  General . 

The  work  presented  here  was  performed  in  three  phases. 
In  the  first  phase  model  velocity  and  drag  were  measured 
experimentally  and  dye-flow  photographs  were  taken  of  the 
model  stern  section  for  a  systematic  series  of  tests.   In 
the  second  phase  theoretical  calculations  were  made  to  de- 
termine the  point  of  separation  on  the  body  which  could  be 
related  to  that  observed  in  the  photographs.   This  phase  in- 
volved the  calculation  of  the  turbulent  boundary  layer  dis- 
placement and  momentum  thickness,  the  pressure-velocity  dis- 
tribution over  the  resulting  corrected  body  and  the  use  of 
these  quantities  in  two  separation  criteria   for  turbulent 
flow.   The  third  phase  was  the  analysis  of  the  results  of 
phases  one  and  two  to  determine  and  relate,  if  possible,  the 
effects  of  the  independent  variation  of  speed,  depth  of  sub- 
mergence, and  stern  shape  on  the  observed  separation  point 
location  on  the  body  as  well  as  to  compare  and  relate  the  ob- 
served point  to  the  calculated  point  of  separation. 

B.  Experimental  Phase. 

In  this  phase  experiments  to  visually  determine  the 
point  of  separation  on  the  model  and  to  measure  the  total 
drag  on  the  model  were  conducted  at  the  Massachusetts  Institute 


of  Technology  Ship  Model  Towing  Tank  under  the  supervision  of 
the  Department  of  Naval  Architecture,  Drag  and  velocity  data 
and  dye-flow  photographs  were  obtained  for  the  model  at  three 
speeds,  at  three  depths  of  submergence  corresponding  to  sur- 
faced, near  surfaced  and  submerged  conditions,  and  with  five 
varying  shaped  tail  sections  which  modified  the  closing  angle 
of  the  model's  stern. 

General  Description  of  Equipment. 

Figure  I  is  a  sketch  of  the  model  showing  the  variation 
in  stern  section  and  the  waterlines  of  the  three  depths  used. 
Figure  II  shows  a  view  of  the  towing  tank  and  towing  carriage, 
while  Figure  III  shows  the  carriage  control  station  and  in- 
strument at  ion . 

The  model,  #4893,  was  constructed  at  the  David  Taylor 
Model  Basin,  C arde rock, Mar yland,aiad  is  a  blunt  nosed  circular 
cylinder  consisting  of  three  main  parts:  a  blunt  solid  ma- 
hogany nose  section,  a  long  hollow  aluminum  parallel  middle- 
body,  which  contains  the  dye  system,  and  five  solid  mahogany 

stern  sections  of  varying  shape  and  length.   To  insure  opera- 
trip 
tion  in  turbulent  flow  a  1/32  inch  diameter/wire  was  placed 

around  the  nose  of  the  model  2.5  inches  from  the  forward  end. 

The  offsets  and  other  characteristics  for  the  model  and  stern 

sections  are  given  in  Appendix  A.   Figure  IV  shows  the  model 

components. 
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Figure  II  -  MIT  Ship  Model  Towing  Tank 


Figure  III  -  Left  to  right:  carriage  control  console, 
digital  voltmeter  and  speed  measuring  equipment,  Sandborn 
recorder. 


Figure  IV  -  Model  and  five  stern  shapes 


Figure  V  -  Internal  view  of  model  parallel  middlebody  showing 
support  plate,  dye  tubes,  dye  manifold  and  brass  dye  ring. 
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The  dye  system  emits  dye  from  1/16  inch  diameter  holes 
spaced  every  1/4  inch  apart  around  a  circular  ring  located  at 
the  after  end  of  the  parallel  middlebody.   This  system  is  shown 
in  Figure  V.   Dye  was  fed  to  the  model  from  a  tank  mounted  on 
the  towing  carriage  through  1/4  inch  plastic  tubing.   A  solu- 
tion of  water  and  potassium  permanganate  was  used  for  dye. 
Each  of  the  five  stern  sections  can  be  attached  to  the  mo- 
del immediately  aft  of  the  dye  ring.   When  the  model  was  in 
motion  dye  released  from  the  dye  ring  would  stream  over  the 
stern  section  permitting  observation  of  the  flow. 

The  model  was  rigidly  mounted  to  the  towing  carriage 
through  a  bracket  which  allowed  operation  at  any  depth  within 
the  limits  of  the  18  inch  towing  street.   Figure  VI  shows 
the  model  mounted  on  the  carriage.   The  towing  carriage  rides 
on  an  overhead  rail  and  is  driven  by  a  variable  speed  elec- 
tric motor  driving  a  soft  rubber  wheel  in  contact  with  a 
similar  side  rail. 

Figure  VI  also  shows  the  drag  measuring  dynamometer 
located  between  the  4:1  ogive  towing  strut,  connected  to  the 
model  internally,  and  the  2j  inch  diameter  hollow-cylinder 
which  was  connected  to  the  mounting  bracket.   The  dynamometer 
measures  the  axial  force  between  the  carriage  and  the  model 
by  means  of  a  differential  transformer. 

A  length  of  32   feet  was  used  for  each  run.   The  length 
of  run  was  limited  to  this  length  due  to  a  portion  of  the 
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Figure  VI  -  Model  mounted  or.  carriage  showing  towing  dynamometer, 
support  bracket,  dye  tubes  and  towing  strut. 


Figure  VII  -  kx5   view  camera  and  strobe  light  arrangement 
at  tank  window. 
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carriage  rail  beyond  this  length  being  badly  distorted  and 
thereby  causing  undue  carriage  vibration.   The  length  of  run 
started  at  the  beach  end  of  the  tank  and  extended  to  the  end 
of  the  first  of  the  10  foot  by  4  foot  glass  windows  located 
in  the  tank's  side  wall.   A  4x5  view  camera  was  set  up  at 
this  window  along  with  two  high  intensity  strobe  lights  which 
were  tripped  by  the  carriage  striking  a  contact  on  the  over- 
head rail  as  it  reached  the  proper  position  in  the  window  near 
the  end  of  the  run.   A  white^  light  reflecting  background  was 
installed  behind  the  model  because  the  distance  to  the  far 
side  of  the  tank  was  too  great  to  achieve  satisfactory  expo- 
sures.  Figure  VII  shows  the  photographic  set-up  that  was 
used.   Each  stern  section  was  marked  at  1/2  inch  intervals 
along  its  length  as  an  aid  in  determining  the  position  of 
separation  from  the  photographs. 

The  carriage  speed  was  set  by  the  manual  speed  control 
dial  on  the  carriage  console  since  the  automatic  speed  con- 
trol was  inoperative.   Speed  data  for  each  run  was  recorded 
on  tape  with  a  Sandbcm  Type  150  recorder  and  type  1300 preampli- 
fier.  Drag  force  was  read  out  on  a  type  405  Hewlett  Packard 
D.C.  digital  voltmeter  with  input  from  the  type  1100  Sandborn 
preamplifier a   The  voltmeter's  purpose  was  to  attenuate  the 
high  frequency  components  of  the  carriage  vibrations  sensed 
by  the  dynamometer  as  drag  fluctuations  and  to  allow  a  volt- 
age proportional  only  to  model  drag  to  be  read  out.  A  more 
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detailed  description  of  the  towing  tank  and  associated  equip- 
ment, model  dye  system  and  dynamometer  is  given  in  Appendix  B. 

Run  Procedure 

Two  runs  were  conducted  for  each  condition  of  stern 
section  depth  and  speed.   On  each  run  a  photograph  was  taken 
of  the  dye-flow  over  the  stern  section  while  speed  and  drag 
were  recorded  simultaneously  on  the  Sandborn  recorder  and 
digital  voltmeter.   The  tank  was  allowed  to  settle  between 
each  run.   Polaroid  pictures  were  taken  occasionally  to  check 
exposure  and  dye-flow  definition.   Only  five  of  the  dye  holes 
in  the  dye  ring  were  used,  two  each  on  the  top  and  the  bottom 
and  one  on  the  axis  of  the  body.   In  turbulent  flow  the  dye 
present  from  the  other  holes  tended  to  obscure  the  separation 
point.   For  laminar  flow  the  use  of  the  entire  ring  was  satis- 
factory.  Figure  VIII  shows  both  turbulent  and  laminar  flow. 
The  unused  holes  were  masked  off.   All  screw  holes,  etc,  on 
the  model  were  filled  smooth  with  a  waterproof  putty  to  insure 
as  little  flow  disturbance  as  possible. 

The  model  speed  was  easily  obtained  from  the  Sanborn 
trace  and  a  speed  calibration  curve;  however,  the  drag  could 
be  obtained  from  the  digital  voltmeter  only  by  averaging  the 
reading  over  the  run.   This  was  done  by  recording  about  ten 
values  of  the  voltmeter  reading  during  the  last  third  of  the 
run  and  averaging.   This  average  reading  was  then  compared 
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to  the  "zero"  drag  reading  taken  prior  to  each  run  to  determine 
a  voltage  difference  reading  proportional  to  drag.   A  calibra- 
tion curve  was  determined  by  applying  a  known  static  load  to 

the  model  and  observing  the  voltmeter  reading.   Using  this  curve 

** 
and  the  voltage  difference  reading  drag  force  was  obtained. 

C.  Theoretical  Phase. 

In  this  phase  a  theoretical  position  of  separation  on 
the  body  was  established  after  correcting  the  body  for  the  ef- 
fects of  boundary  layer  growth.   Throughout  this  phase  it  was 
assumed  that  the  flow  was  turbulent  from  the  nose  of  the  body 
and  that  the  1/7  power  law  holds  for  the  velocity  distribution 
within  the  turbulent  boundary  layers. 

Calculation  of  Turbulent  Boundary  Layer  Characteristics. 

As  indicated  in  the  general  discussion  of  the  procedure 
the  value  of  the  displacement  thickness  on  the  body  and  mo- 
mentum thickness  on  the  body  are  necessary.  Eckert  [4]  has 
evaluated  the  effects  of  transverse  curvature  on  the  turbu- 
lent boundary  layer  on  a  circular  cylinder  in  compressible 
flow.   These  effects  are  the  same  as  reported  by  Yu  \s] 
and  Landweber  ((f) ,  which  are : 


**  The  drag  measurements  were  unsatisfactory.   Appendix  E 

has  a  more  detailed  discussion  relating  to  these  measurements. 
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(1)  Boundary  layer  thickness  and  displacement  thickness 
are  less  than  on  the  flat  plate. 

(2)  Momentum  thickness  and  local  and  mean  skin  friction 
coefficients  are  higher  than  flat  plate  values.  _ 

(3)  The  effects  of  (1)  occur  when  the  ratio  of  o/r-^O.l. 

(4)  The  friction  coefficients  deviate  from  flat  plate 
values  only  when  ratio  £/y  &   1.0. 

Eckert  develops  simple  equations  which  relate  the  value 

of  6  )    o  ,  and  ©  on  the  body  to  values  of  the  same  parameters 

on  a  flat  plate.   The  assumption  is  made  that  the  velocity 

profile  within  the  boundary  layer  conforms  to  the  l/7th  power 

law  distribution  of  U  - (Y)  and  that  Blasius1   Law  relating 

UooW 
wall  shear  stress  as  inversly  proportional  to  the  fourth  root 

of  the  boundary  layer  thickness  are  uneffected  by  the  curva- 
ture.  The  assumption  is  also  made  that  the  boundary  layer  is 
turbulent  from  the  beginning  and  it  starts  from  the  nose  of 
the  body  with  zero  thickness. 

The  relationships  developed  for  incompressible  flow  are 
given  below  and  the  reader  is  referred  to  [4]  for  the  deriva- 
tion details  which  are  straightforward. 
Relationship  for  boundary  layer  thickness 

6   -  A A4'5    (i) 

where 

6  pi  -  0.38  X       (ia) 
(RxT772 

Although  (1)  is  implicit  in  o  it  can  be  easily  solved  by  using 
£  pi  as  a  first  try  and  iterating  the  results. 
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Relationships   for     Displacement   thickness. 


>  *.  i  +<S    --1/1  +  2<?  (l- A„V   i2fl-A       ^1  (2) 

iA,pl^/3r)4/5 

where,  *  *  r 

A      „i    =  O-   P1  "     1  (2a)     £  m"  2£l  <2b> 

^2npl   a^TSmU       2WT-1        (2c)  5t    =a- 

Relationships    for  Momentum  thickness. 

where , 

r>  C     =    7^      -     221       "  ^  0  pl  (3a) 

P1  ?i+l         ?^+2  ^+l(rt+2)      3  Pl 

-ft    2-n,„l     -J2~      ,         (77.»2)        -  Vis  (3b) 


r   * 
The  value  of  ^>   calculated  from  (2)  was  applied  as  a 

correction  to  the  body  surface  after  the  manner  described  in 

(7}  as  is  shown  in  Figure  X  taken  from  [7j . 
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J3o£>y    /QX/S 


Figure  X  -  Equivalent  Body 
Profile  with  Displacement  Thickness 


Thus  the  corrected  body  will  have  a  radius  given  by 
V  corr  -  r  +  Cos  t     &  (4) 

For  use  in  the  calculation  of  the  pressure  velocity  dis- 
tribution discussed  in  the  next  section  this  radius  was  norma- 
lized as  follows: 

V   corr  -  I^corr (4a) 

norm   body  length/2 

An  IBM  7090  computer  program,  Program  II,  was  written  to 
evaluate  equations  (1,2,3,4  &  4a)  at  16  points  on  each  of  the  5 
bodies  for  9  speeds.   These  16  points  correspond  to  the  loca- 
tion of  the  16  ordinates  required  for  the  Gaussian  16  point 
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quadrature  rule  use  in  the  determination  of  the  pressure  veloci- 
ty profile.   The  values  of  body  radius  and  the  angle  o  were 
picked  off  at  each  of  the  16  points  from  full  scale  drawings  of 
the  5  bodies  prepared  by  the  authors.   A  tabulation  of  this 
data  is  contained  in  Appendix  A.   Details  of  the  computer  pro- 
gram are  in  Appendix  C, 

Calculation  of  Pressure  Velocity  Distribution 

Potential  Flow  Pressure  and  Velocity  Distribution  was 
calculated  according  to  a  method  developed  by  Landweber  [&]  . 
Two  types  of  solutions  are  given  for  obtaining  the  steady,  ir- 
rotational  ,  axisymmetric  flow  of  an  inviscid,  incompressible 
fluid  about  a  body  of  revolution.   The  first  method  involves  an 
iterative  type  solution  for  Fredholm  integral  equation  of  the 
first  kind  and  is  not  used  in  this  work.   The  second  method  ob- 
tains  the  velocity  directly  as  the  exact  solution  of  an  inte- 
gral equation  of  the  first  kind  by  an  arithmetic  procedure. 
Using  the  second  method  as  outlined  in  [&]   the  basic  problem 
is  to  solve  for  the  velocity  potential  0  which  is  dependent  on 
x,  y  and  satisfies  Laplace's  equation  in  cylindrical  coordinates 

JL  fyggV-a    (v  MS  -  o  <*> 

ox  v   ox/      c)  y  v.       ay/ 

where  the  equation  of  the  normalized  body  profile  is  in 

the  form:  y  -  f  (t0)   (Figure  XI  shows  method  of  normalizing  the 

body) 

and  u  -  -  a.%        ,  as   -  -  o  0 

c)x  dy 
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Figure  XI  -  Normalized  Body  Profile 


There  also  exists  a  Stokes  stream  function  ^(x,y) 
which  satisfied  the  equation: 

^2£  +  g>2^,  J.  %p  (6) 

^x2      3y2      y  ^y~ 

For  uniform  flow  of  velocity  IX    parallel  to  the  x-axis 

0  -  -Z/x  and^-  -  JLZy2  at  infinity 

The  following  boundary  condition  on  the  body  must  also  be 

satisfied: 

(a0\-   -  V  cos/3 

Where  V  »  velocity  of  the  body  and  £»-  the  angle  between  the 

outward  normal  to  the  body  and  the  x-axis. 

Derived  from  Green's  Theorem  the  equation: 


//< 


"n-'s  -      -  a£-  « 


(7) 
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where  0  and  w  are  any  functions  harmonic  in  the  region  exterior 

to  a  given  body  and  vanishing  at  infinity,  is  used  to  obtain 

an  integral  equation  for  axisymmetric  flow  about  a  body  of 

revolution. 

Substitution  of:   dS  -  2  "77"  y  ds  and 

d0  -  -  sin  X         into     (7)   gives: 
dn 

I     y0  dw  ds  -  -  I  y  w  sin  )f  ds        (8) 
Jo  3n       JQ 

y  -  ordinate  of  meridian  section  of  the  body  which  is 

moving  in  the  negative  x  -  direction  with  unit  velocity. 

ds*-  element  of  arc  length  along  body 

6   ■  angle  of  the  tangent  to  the  body  with  x  -  axis 

2P-  perimeter  of  a  meridian  section 

With  W  as  an  axisymmetric  potential  function  and 

T  (x,y)  the  corresponding  stream  function 


y  dw  -  6J£ 
dn   ds 


and   _P 

y0  dw  ds  -  0(j/ 

<o  dn 


P   sP 

-  J  V   ^0  ds 

0  Jo       ds 


If  LC   is  total  velocity  along  a  body  when  a  stream  of 

unit  velocity  is  superposed  in  the  positive  x  -  direction 

H     ---dp;  +  costf 
ds 

Since  dx  -  ds  cos^  and  dy  -  ds  sin  &  equation  (8)  can 

be  written  as: 
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p     /P 

Upds*l     (fdx   -  yw  dy)  -0ff 
■lo 


P 

o 


(9) 


but   y  -  3w  -  -  d£     since  w  and  ^  are  corresponding  potential 

2T?.     Jy 
and  stream  functions  and  ^  dx  -  ywdy   is  an  exact  differen- 
tial defining  a  function  a2  (x,y)  such  that 


9£ 


ay 


From  this  and  y  bw  -  g)£     we  obtain 

a  x2     3  y2      y  "ST 

which  is  identical  to  (6)  and  is  the  equation  satisfied  by 
Stokes  stream  function.   It  can  also  be  verified  that  w  and  P 
are  corresponding  axisymmetric  potential  and  stream  functions 


which  satisfy  the  equation:  9<^  -  -  y  0<# 
Writing  (9)  in  terms  of  _d 

q-q  dSi  -  (n  -  0  sji  \ 


a* 


£>y  / 


2f?-  y   £>0 
<2y  dx 


(10) 


o 


If  JT  is  the  stream  function  of  a  point  source  of  unit 
strength  situated  at  x  ■  t  •  the  x  -  axis 


-ft 1  +  x-t 

Then  c)ft 


2 


"Sx 

('vSl-  0  ^ft   ") 


,  Jl   -  [(x-t)2  +  y2J     (11) 

,  and  if  y  vanishes  at  both  limits, 
P 


df 
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(11)  now  becomes: 

U   (x)y2  (x)c/s    -  l  (12) 

2  Jt  3 
o 

which  is  an  integral  equation  of  the  first  kind  in  which  the 

unknown  function  is  LL  (x)  and  the  kernel  is  y  /oj?*      When  XZ 

is  determined  the  pressure  distribution  is  given  from 


Bernoulli's  equation. 


2  A  .^  -1 


E_  -  1   -  (u2  +  ^2)  »  1-U 

q 

where  q  is  the  stagnation  pressure. 

Landweber  uses  an  iteration  procedure  to  solve  (12) 
as  follows: 

If  the  iterations  are  expressed  in  terms  of  error  functions 
En  ^^  at  the  Point  x  -  t  ,  first  obtain 
E1  (t)  from    En  (t)  -  1  -  J  U n  (x)  y2  (x)  ds 

by  using  the  first  approximation  [J.  (t)  -  (l+k)  cos  a   (x) 
where  \a±   -  longitudinal  virtual  mass  coefficient.   Then  suc- 
cessive iteration  E2,  E3  E   are  obtained  from 

Vl  (t)  "  En(t)  "  *f  En  (x)  y2(x)   dx 

where^0    ,  ■•/•       are  nose   and  tail    abscissaeof   the  body   and 

finally  |]n+|(t)    is   obtained  from 

n 
Un+i    (t>   "    LLt  <t)  +  cos     %    (t)      21      *i    Ct)  (13) 
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To  evaluate  integrals  of  the  form 

x/  2 

E    (x)   y      (x) 

x0  St  3 

use     k      (x,t)   *  -f  (x)  (14) 

Q(x-t)2  4  f  (x)  J^2 

k'    <V>   "  g    (  y-i  t)  (15) 

[( *-t)5  +  Jc^Tpi 

where  y  -  j  00  is  the  equation  of  the  given  body  profile 

and  y2  -  g(x^  t))  is  the  equation  of  the  prolate  spheroid 

whose  ends  coincide  with  the  ends  of  the  given  body,  and 

which  intersects  the  body  at  x  -  t,  i.e.} 

g(*,  t)  -/(t)   (x-ypHx!  -x)  (16) 

(t-  x0)CXi  -  t) 

For  the  spheroid  the  length  -  diameter  ratio  71  is  given 

by      X2  -  (t  -  xp)  (x,-  t)  (17) 

and  from  this  the  virtual  mass  coefficient  K4  (t)  can  be  ob- 
tained from: 


k4   -   Xln(X-4/>T^l)  -  jrf-   1 
1      P  fjS   _f  1  >lSCM  2  -1)  <18> 

E±(t)is  then  obtained  from 

B±    (t)  -  1  -  1+kJ  JK  (x^t)   -k'   (x>t)]   dx  -     1  +  ki     (19) 

Z,       J  J.  HI  |C«  kX  / 

and  finally 


En+1    (t)  "  E    (t)-Hci  E 
I  +  k 


n(t)   -J      /    k(X,  t)JEn(x)-En(t)]dx 
A<o  <20> 
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In  using  Landweber*s  method  the  velocity  is  obtained 
at  16  points  along  a  body,  from  t  -  -1  to  t  -  +  1,  correspond- 
ing to  the  Gaussian  values  (<^)  for  the  16-point  guadrature 
rule.   In  applying  the  rule  the  integrands  in  (19)  and  (20) 
are  evaluated  at  the  16  Gaussian  abscissae  Xj  -  £  .  for  each 
of  the  16  values  of  t±-  %  ±    .   Gaussian  weighting  functions 
are  used  to  obtain  K  .  .  =  R  j  k  (X^  ti)  and  Kj  ^  -  Rjk 
(Xj  t   t±)   which  are  used  in  evaluating  (19)  and  (20).TIn(t)  is 
then  given  by  (13)  and  ?/&   by  Bernoulli's  equation  in  the  form 

P/q  -  1  -Un2  • 

In  the  application  of  Landweber*s  method  of  determining 

the  pressure-velocity  distribution  on  a  body  of  revolution  to 
the  models  used  in  this  thesis  the  corrected  body  profile  as 
calculated  in  the  previous  section  was  used.   The  function 
Y*  -jf (x)  for  the  normalized  corrected  body  profile  was  not 
determined  because  even  if  it  is  determined  the  end  result  is 
still  a  point- wise  distribution  of  ^/Lk>  .   Instead  of  determin- 
ing the  polynomial  describing  the  body  profile  the  values  of 
the  normalized  corrected  body  radius  at  each  Gaussian  ordinant 
computed  in  Program  II  were  used.   The  angle,  #  ,  was  again  de- 
termined from  full  scale  drawings  of  the  model  and  stern  sections 

An  IBM  7090  computer  program,  Program  III,  was  written  to 
determine  the  pressure-velocity  distribution  on  the  body  cor- 
rected for  boundary  layer  effects  using  Landweber*s  method.  The 
program  calculates  Cp,  L£/Xlo©>  absolute  velocity  and  gage  pres- 
sure on  the  body  at  16  points  on  each  of  5  bodies  for  3  different 

24 


speeds  at  each  of  3  different  depths  of  submergence.   The  9 
speeds  used  are  the  same  as  those  used  in  Program  II  and  are 
the  average  of  the  speeds  listed  for  each  run  and  its  dupli- 
cate, i.e.,  Run  1  and  Run  1A,  in  the  tabulation  of  experimen- 
tal data  in  Appendix  D.   Thus  the  outputs  of  the  two  programs 
are  a  compatible  package  for  application  in  the  separation 
criterion  to  be  evaluated.   This  output   is  in  turn  compati- 
ble with  the  dye-flow  photographs  since  the  same  conditions 
apply*   A  test  solution  for  Program  III  was  run  on  the  ex- 
ample worked  out  in  [7] . 

Separation  Criteria  . 

Stratford  [2]    has  developed  a  rapid  method  for  the  pre- 
diction of  flow  separation  in  turbulent  flow  from  an  approxi- 
mate solution  of  the  equations  of  motion.   The  equations  are 
integrated  by  a  modified  "inner  and  outer  layer  solution" 
technique  developed  by  Stratford  in  1954  for  laminar  flow 
boundary  layers.   The  final  solution  is  given  by  equations 
(21)  through  (24)  which  relate  X,  Cp,  and  dCp/dX  at  the  sepa- 
ration position  and  may  be  applied  to  any  given  pressure  dis- 
tribution.  The  algebra  and  procedure  for  obtaining  the  sepa- 
ration condition  are  clearly  detailed  in  [2]  and  only  the  re- 
sults will  be  given  here. 

The   resulting  implicit  equation  for  the  separation  po- 
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sition,  X,  is  -  _  1 

(2Cp)i<°-2>  (X  _g)  -  1.06^  (10-6RS)  ^   (21> 

for  Cp  -  — - 

n+1 

Where:  R  is  the  Reynolds  number  based  on  value  of  distance 

X  and  freest ream  velocity,  Uco. 

n  =  log  10  Rs  (22) 

and^-  0.66   for  dja  <  0  (23) 

dX2 

j9    -  0.73   for  d2p   g  °  <24> 

dX^ 

2 

with  the  resultant  value  of  dp  being  that  immediately  prior 

dX2 

to  separation,   This  formulation  has  a  range  of  uncertainty 

of  10%  due  to  its  inability  to  fully  account  for  d  p/dX  in 

the  case  of  turbulent  boundary  layers.   However,  the  rapidity 

of  the  method  makes  this  range  of  uncertainly  acceptable  in 

this  case,   Granville  in  [3]  has  brought  forward  another  even 

simpler  separation  criterion   It  was  apparently  originally 

developed  by  Bu^i  in  1931  at  Zurich,   The  criterion  states 

that  separation  will  occur  if  the  following  relation  holds: 

9  Uoo  v  du y  L   <   -0,006  (25) 

ux~  X  dX*U3o"  - 

Where;  8 =  momentum    thickness  on  the  body  at  separation 
point  (Calculated  using  equation  C3^>  )  and  L  -  body  length. 

since  L  and  Uoo  are  constants  for  each  body  at  a  parti- 
cular speed 

/  U\ 
du  =   d  — —    1T 


(i) 
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Equation  (25)  can  thus  be  written 

•*s*2B.  x<»5r)--0-006  (26) 


U 


Wl 


Unfortunately  the  information  regarding  this  criterion  is 
limited  to  that  obtained  by  the  writers  from  Course  13.012  in 
the  Department  of  Naval  Architecture  at  M.I.T.  and  the  meager 
yet  reliable  reference  to  it  in  [3]  .   A  range  of  uncer- 
tainty for  the  criterion  is  not  known.   It  is  intended  to  be 
used  as  a  criterion  to  be  used  in  scaling  separation  from 
model  to  ship  scale. 

With  the  velocity  profiles,  pressure  coefficients,  and 
momentum  thickness  values  on  the  body  determined  as  outlined 
in  the  previous  sections  it  is  possible  to  evaluate  both  of 
these  criteria  and  see  how  they  will  compare  with  the  experi- 
mentally observed  points  of  separation.   This  was  done  by 
hand  calculation  for  a  few  cases.   Sample  calculations  are  given 
in  Appendix  F.   The  results  of  such  calcuations  are  in  the 
RESULTS  section  of  this  thesis, 

D.  Analysis  Phase. 

The  data  obtained  from  the  experimental  phase  relating  to 
drag  was  analyzed  by  utilizing  an  IBM  7090  computer  program, 
Program  I,  to  reduce  the  data  to  coefficient  form.   The  details 
of  this  computer  program  are  given  in  Appendix  C.   The  analysis 
of  the  experimental  drag  data  contained  in  Table  D-I  of  Appen- 
dix D  showed  the  data  to  be  unreliable  and  faulty.   Appendix  E 
contains  some  of  the  results  obtained  from  the  computer  program 
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and  a  discussion  relating  to  these  results.   The  second  step 
in  analyzing  the  experimental  data  was  to  study  the  dye-flow 
photographs  and  determine  from  them  the  point  of  separation 
and  the  body  angles  with  respect  to  the  body  axis  at  the  point 
of  separation.   The  point  where  the  dye-flow  appeared  to  a- 
bruptly  depart  from  the  body  profile  was  taken  to  be  that  of 
the  point  of  separation.   Only  the  flow  along  the  keel  of  the 
model  was  used  in  determining  the  location  of  the  point  of 
separation  because  it  was  felt  that  the  presence  of  the  towing 
strut  forward  would  influence  the  flow  over  the  top  of  the 
model.   The  authors  each  made  independent  determinations  of 
the  separation  position  on  each  run.   Differences  which  arose 
in  interpreting  the  location  of  the  separation  point  were 
then  compromised.   The  determinations  were  also  made  prior  to 
the  calculation  of  the  separation  position  from  the  two  sepa- 
ration criteria  so  that  the  location  of  these  theoretical 
points  would  not  influence  the  interpretation  of  the  photographs 
In  this  manner  it  was  felt  that  as  high  as  possible  degree  of 
precision  would  be  obtained  in  interpreting  the  photographs. 
There  was,  as  was  to  be  expected,  considerable  judgment  in- 
volved in  establishing  the  separation  point  in  several  of  the 
pictures.   A  sampling  of  typical  photographs  from  which  the 
separation  point  was  determined  is  shown  in  Figure  XII.   The 
location  of  the  observed  separation  point  for  each  run  is  con- 
tained in  Table  D2  in  Appendix  D.   The  angle  between  the  body 
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Run  1 


Run  2 


Run  3 


• — Hr-.i  1 1 1 1 1 y  n* 


Run  4 


Run  5A 


Run  6 


-- 


Run  7  Run  8  Run  9A 

Ao  Sequence  of  Runs  on  Body  3e 


Run  13A 
Fiody  1 


Run  12A 

Body  1 


Run  20A 
Body  2 


Run  35 A 
Body  U 


Run  29 
Body  4 


Run  38A 
Body  5 


B.  Runs  on  Cther  Bodies. 
Figure  XII  -  Typical  Dye-flow  Photographs 
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surface  and  the  body  axis  at  the  point  of  separation  was  de- 
termined graphically  by  plotting  the  separation  point  loca- 
tion on  full  scale  drawings  of  the  bodies.   The  movements  of 
the  point  of  separation  were  analyzed  to  see  what,  if  any, 
correlation  could  be  made  between  these  movements  and  changes 
of  speed,  depth  of  submergence  or  body  shape.   The  angle 
at  which  separation  occurred  was  similarly  analyzed. 

The  third  step  in  the  analysis  was  to  calculate  the 
point  of  separation  by  each  of  the  two  separation  criteria 
using  the  results  of  computer  programs  I  and  II.   Samples 
of  the  output  from  these  programs  are  contained  in  Appen- 
dix G.   Sample  calculations  for  each  of  the  criteria  are 
given  in  Appendix  F.   During  these  calculations  the  criti- 
cality  of  the  various  parameters  in  the  criteria  were 
noted.   The  results  of  these  calculations  were  compared 
to  the  observed  point  of  separation  noting  location  and 
angle  of  body  surface  at  the  point  of  predicted  separation. 
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Ill  RESULTS 

The  results  of  the  procedures  described  in  the  pre- 
vious sections  are  presented  in  this  section.   These  results 
are  discussed  and  conclusions  are  drawn  from  them  in  the 
next  section  of  the  thesis. 

The  results  are  presented  primarily  in  graphical 
form.   Figures  XIII  through  XVII  show  the  relationship  be- 
tween the  pressure- velocity  distribution  over  the  body  and  the 
location  of  the  observed  separation  points.   Figures  XVIII 
through  XXII  show  the  location  of  the  experimental  and  calcu- 
lated points  of  separation.   The  body  profiles  corrected  for 
boundary  layer  effects  are  also  shown  here.   Table  I  con- 
tains a  breakdown  of  the  experimentally  observed  separation 
points  on  each  body  by  depth  of  submergence,  and  speed.   The 
body  angle   at  separation  is  also  shown. 
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TABLE  I 
Results  of  Studies  of  Dye  Flow  Photographs  and  Separation  Criteria 


BODY  #1 
Code    Angle 


\ 
B 
C 
E 
F 


-  240 

-  32° 

-  420 

-  6° 

-  47° 

BODY  #2 
Angle 


B 
C 
D 
E 
F 


- 25^ 

-  30° 

-  36° 

-  42° 


,0 


-  13' 

-  38° 
BODY  #3 

Angle 


XT 


B 
C 

E 
F 


B 
C 
E 
F 


-  23 

-  28° 

-  33° 

-  10° 

-  46° 
BODY  #4 

Angle 

-  2fio 

-  38° 

-  49° 

-  6° 

-  36° 
BODY  #5 

Angle 


B 
C 

E 
F 


-  23^ 

-  30° 

-  60° 

-  5° 

-  23° 


.955" 

.970 

,978 

.94o 

,982 


Points  (Speed-Depth) 
5  @~LM,  2  %   LD,  2  @  MD,  2 
3  ©  MM,  2  (&   HM,  1  ^  MS 
2  @  HD,  2  41  HS,  2  ©  MS 
BURI  CRITERION 
STRATFORD  CRITERION 


LS 


x/l 

Points  (Speed-Depth) 

.967 

1  m   LD,  1  m  LM 

.976 

2  @  MM,   1  43  HM,  1  g  LD, 

.983 

2  t  MD,  1  @   HD,  1  41  MS,  .1 

,989 

2  m   HS,  1  @  MS 

.930 

BURI   CRITERION 

.986 

STRATFORD  CRITERION 

X/L 

Points  (Speed  -  Depth) 

.967 

2  m  ms,  2  m  ls,  1  43  hs 

.975 

2  @  LM,  2  43  MM,  1  43  HM,  1 

.982 

2  «  HD,  2  €  MD,  2  €  LD   1 

.905 

BURI   CRITERION 

,992 

STRATFORD  CRITERION 

2   @  LS 
@  LM 


42  HS 
@  HM 


Points    (Speed 
4l~MD 


Depth 


2  42  LD, 
2  43  LS,  2  €  MS,  1  m   HD9 
2  €  HM,  2  €  MM.  2  €  HS, 
BURI   CRITERION 
STRATFORD  CRITERION 


m  md 
m  lm 


1  €  LM 


Points  (Speed-Depth) 
All  speeds  at  deep  depth 
All  speeds  at  medium  depth 
All  speeds  at  surface 
BURI   CRITERION 
STRATFORD  CRITERION 


Separation  points  are  plotted  on  body  profiles  in 
figures  XVIII  through  XXII. 

Point  designations  are  the  same  as  outlined  in  Table  D-II 
of  Appendix  D;  e,g.   2  43  HM  means  two  points  were  plotted  at 
that  location  and  angle  on  the  body  corresponding  to  two  runs 
at  the  high  speed  -  medium  depth  condition. 
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IV  DISCUSSION  OF  RESULTS 

Accuracy  of  Observations 

The  accuracy  achieved  in  determining  the  location  of 
the  separation  points  from  the  dye-flow  photographs  was 
limited  by  the  ability  of  the  authors  to  pick  off  the  point 
where  the  flow  abruptly  departed  from  the  body  surface  and 
in  their  ability  to  interpolate  between  the  J  inch  markings 
on  the  bodies.   In  some  cases  the  point  of  separation  was 
partially  obscured  by  excessive  amounts  of  dye  while  in 
other  cases,  especially  at  the  higher  speeds,  a  portion  of 
the  dye  appeared  to  be  carried  beyond  the  separation  point 
by  flow  outside  the  boundary  layer.   In  the  surface  condition 
there  was  some  evidence  of  dye  being  drawn  up  into  the  area 
aft  of  the  body  toward  the  free  surface.  The  accuracy  in 
observing  the  point  of  separation  in  turbulent  flow  can  not 
be  equal  to  that  obtainable  for  laminar  flow,  but  improvement 
in  this  accuracy  could  be  gained  by  using  a  larger  model  and/or 
using  a  more  sophisticated  means  of  photographing  the  flow 
such  as  telephoto  lens  or  underwater  motion  pictures  to  ob- 
tain greater  definition  of  the  flow  about  the  body. 

Effect  of  Body  Closing  Angle 

In  examining  the  positions  of  the  observed  points  of 
separation  on  the  various  bodies  as  shown  in  Figures  XVIII 
through  XXII  it  was  noted  that  the  points  occurred  in  groups 
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located  at  three  or  four  positions  on  the  bodies.   Table  I 
shows  the  details  for  these  groupings.   The  center  of  each 
group  was  used  as  the  measuring  point  of  the  body  angles.   The 
degree  of  accuracy  obtained  in  measuring  these  angles  is  a 
function  of  the  accuracy  achieved  in  the  initial  location  of 
the  point  of  separation  already  discussed  and  the  errors  en- 
voi ved  in  measuring  the  angle  from  the  full  scale  drawings. 
The  two  factors  result  in  a  possible  error  of  -  3  degrees. 

The  initial  group  of  separation  points  on  all  five 
bodies  tested  occurs  at  a  body  angle  of  between  23  and  28 
degrees,  averaging  about  25  degrees.   There  is  no  apparent 
orderly  correlation  among  the  five  bodies  to  account  for  this 
variation  except  that  it  is  due  to  experimental  error.   This 
results  in  the  conclusion  that  the  initial  point  of  separa- 
tion occurs  at  approximately  the  same  angle  for  all  bodies 
regardless  of  the  amount  of  streamlining.   A  critical  angle 
of  25  degrees  thus  appears  to  be  established  for  initiation  of 
separation.   The  rate  at  which  the  body  angle  is  changing  does 
not  appear  to  effect  this  critical  angle  because  bodies  1  and  5 
which  have  vastly  different  rates  of  change  of  body  angle  have 
nearly  the  same  critical  angle.   A  systematic  check  over  a 
greater  range  of  body  shapes,  speeds  and  depths  of  submergence 
is  warranted  to  check  this  conclusion. 

The  maximum  value  of  the  angles  in  the  region  over 
which  the  separation  points  occurred  varied  from  about  42  degrees 
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to  60  degrees  showing  that  the  observed  point  of  separation 
did  fluctuate.   At  these  higher  angles,  however,  the  speeds 
were  generally  toward  the  high  end  of  the  speed  range  and  the 
fluctuations  may  be  caused  by  dye  being  carried  aft  by  flow 
outside  the  boundary  layer,  especially  on  body  5  where  the 
point  of  separation  was  most  difficult  to  determine.   Be- 
cause of  this  possible  effect  these  angles  are  not  con- 
sidered as  significant  as  the  initial  separation  angles. 
The  above  comment,  notwithstanding,  these  higher  angles  do 
point  up  that  under  certain  circumstances  separation  of 
flow  may  not  occur  until  the  critical  angle  has  been  ex- 
ceeded by  a  considerable  amount. 

■ 

Evaluation  of  Separation  Criteria 

Buri's  theoretical  separation  criterion  equation  (26), 
consistently  predicts  the  separation  point  well  forward  of 
both  the  experimentally  observed  points  and  the  point  of 
separation  established  by  Stratford's  criterion,  equation 
(21).   For  Buri's  criterion  the  body  angle  values  all  fall 
between  5  and  13  degrees  with  no  apparent  consistency  of 
variation,  except  that  when  body  1  is  eliminated  the  angles 
decrease  from  body  2  to  body  5,  thus  indicating  that  the 
separation  point  moves  forward  with  increase  in  streamlining. 
This  is  inconsistent  with  the  experimental  results  and  the 
general  concepts  relating  to  separation  which  maintain  that 
streamlining  should  decrease  the  extent  of  the  separation  zone, 
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not  increase  it.   This  criterion  is  very  sensitive  to  slope  of 
the  curve  of  'Uoo  versus  'L  as  plotted  in  Figure  F-I,  Appen- 
dix F,  and  relatively  insensitive  to  model  speed  that  enters 
through  0  in  (26).   The  criterion  predicts  separation  to  occur 
ij.ust  as  the  value  d\    Uoo) /d  \      L)    becomes  negative,   A  slight 
error  in  plotting  this  curve  or  in  using  it  to  solve  (26) 
could  result  in  a  larger  error  in  the  separation  point;  how- 
ever, the  magnitude  of  the  possible  error,  -0.2  inch,  is  not 
sufficient  to  improve  the  agreement  between  Buri  and  the  ob- 
served points. 

On  the  other  hand,  Stratford*s  criterion  gives  much 
closer  agreement  with  the  experimentally  observed  points  of 
separation,  as  it  consistently  predicts  separation  either  with- 
in or  slightly  aft  of  the  observed  region  of  separation.   The 
values  of  the  body  angles  at  the  predicted  separation  points 
vary  from  23  to  47  degrees.   Again  there  is  no  apparent  con- 
sistency of  this  variation.   Depending  on  the  shape  of  the  curve 
of  Cp  versus  X,  as  plotted  in  Figure  F-II,  Appendix  F,  this 

criterion  is  sensitive  both  to  the  value  of  Cp  and  the  value 

d  'C  \ 
of  I   p'/dx  and  also  insensitive  to  model  speed  that  enters 

through  R_  and  n  in  (21),   Where  the  slope  is  steep  the  value 
of  Cp  is  the  controlling  factor  while  for  a  shallower  slope 
both  Cp  and  d(CpVdx  control.   The  criterion  essentially  pre- 
dicts separation  to  occur  very  shortly  after  Cp  becomes  posi- 
tive.  Again  errors  in  plotting  or  interpreting  Figure  F-II 
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would  not  effect  the  predicted  point  of  separation  by  more  than 
±.0.1  inch. 

It  can  be  concluded  that  Buri's  criterion  did  not  give  a 
satisfactory  comparison  to  the  observed  point  of  separation 
while  Stratford's  agreed  much  better  although  it  generally 
indicated  that  the  body  angle  at  separation  was  of  a  larger 
value  than  that  observed.   It  is  to  be  kept  in  mind,  however, 
that  these  criteria  were  analyzed  on  the  basis  of  the  calcu- 
lated profiles  in  Figures  XIII  through  aVII  which  were  de- 
rived for  modified  potential  flow.   It  is  therefore  recom- 
mended that  an  experimental  determination  of  velocity-pres- 
sure profiles  be  made  on  the  bodies  and  that  these  be  applied 
to  the  separation  criteria  before  a  final  judgment  is  madet 

Effect  of  Proximity^  of  Free  Surface 

Table  I  yields  no  apparent  consistent  variation  of  the 
experimental  points  of  separation  on  the  various  bodies  with 
a  change  in  operating  depth. 

Computer  .t'rogram  III  calculates  gage  pressure  or  dy- 
namic pressure  plus  hydrostatic  pressure  due  to  depth  of  sub- 
mergence on  each  body.  The  variation  in  pressure  at  the  keel 
is  from  about  19  lbs.  per  square  foot  at  the  surfaced  condi- 
tion to  about  93  lbs.  per  square  foot  at  the  deep  depth  which 
should  1c  suf/.'icloni;  tc  show  the  effect  of  change  in  depth  on 
the  observed  location  ox  i\ie  separation  point  if  in  fact  such 
an  effect  on  separation  does  exist.   Loth  of  the  separation 
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criteria  are  completely  independent  of  depth  since  it  does  not 
enter  into  any  of  the  variables  contained  in  equations  (21)  or 
(26).   It  is  concluded  from  the  experimental  results  that  hydro- 
static pressure  over  the  range  investigated  has  no  recognized 
consistent  effect  on  the  location  of  the  separation  point. 

It  is  recommended  that  further  experimentation  be 
conducted  with  larger  variations  in  depth.   Perhaps  the  use 
of  a  water  tunnel  or  propeller  tunnel  would  be  better  than 
a  towing  tank  for  this  type  of  experiments  since  a  large 
static  pressure  could  be  applied  to  the  model  in  the  test 
section.   More  typical  ship  shapes  rather  than  body  of  re- 
volution forms  should  be  used  in  conducting  dye-flow  experi- 
ments in  the  air-water  interface  to  show  the  effects  of 
draft  on  separation  over  a  given  form. 

Effect  of  Speed 

Bodies  I  through  4  show  values  of  the  body  angle  at 

the  separation  point  between  23  and  28  degrees  for  the  low 
speed  range,  between  30  and  37  degrees  for  the  medium  speeds 
and  33  to  60  degrees  for  the  high  speeds.   Because  of  the  dif- 
ficulty in  locating  the  separation  point  for  body  5,  no  rela- 
tion between  speed  and  separation  point  could  be  seen  for  that 
body.   The  results  for  bodies  1  through  4,  however,  show  that 
the  separation  point  appears  to  move  aft  on  the  body  as  speed 
increases.   The  reason  for  this  effect  could  be  due  to  the 
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fact  that  dye  was  being  carried  back  along  the  body  outside 
the  boundary  layer  and  misinterpretation  of  the  photographs 

resulted.   At  the  speeds  used  the  Reynolds'  Numbers  were  from 

5 
3.5  to  7  x  10  .   This  is  within  the  doubtful  transition  region 

between  laminar  and  turbulent  flow  and  the  possibility  exists 
that  the  lower  speeds  were  at  the  laminar  end  of  the  transi- 
tion, although  an  examination  of  the  photos  in  Figure  XII 
indicates  a  typical  turbulent  flow  pattern.   As  speed  in- 
creased the  mere  fully  turbulent  flow  region  was  approached 
with  the  added  increased  effect  of  the  trip  wire  with  speed 
thus  causing  the  point  of  separation  to  move  aft  on  the  body, 
as  it  does  in  a  typical  transition  from  laminar  to  turbulent 
flow.   The  two  theoretical  separation  criteria  are  insensi- 
tive to  changes  in  speed  in  the  range  investigated  as  has 
been  previously  discussed. 

The  recommendation  here  is  obvious;  the  experiments 
must  be  conducted  at  higher  speeds  to  further  insure  turbu- 
lent flow.   In  the  present  experiments  the  speeds  obtainable 
were  directly  related  to  the  ability  to  observe  the  dye  flow, 
the  dye  being  dissipated  with  increased  speed  and  almost 
wholly  indistinguishable  beyond  1.75  knots.   The  speed  limi- 
tation was  due  in  part  also  to  large  carriage  vibrations  at 
higher  speeds  than  those  used.   With  a  more  refined  method  of 
turbulence  stimulation  and  dye<=flow  observation  techniques 
higher  speeds  could  be  attained  leading  to  perhaps  more  con- 


49 


elusive  and  plausible  results  relating  the  effect  of  changing 
speed  on  the  separation  point  than  those  herein  achieved. 

Corrected  Body  Profiles  and  Pressure-Velocity  Profiles 

Equivalent  or  corrected  body  profiles  shown  in  Figures 
XVII  through  XXII  are  the  results  of  Computer  Program  II. 
From  the  sample  of  these  results  given  in  Table  G-I  in 
Appendix  G  it  can  be  noted  that  there  is  only  a  very  slight 
variation  in  corrected  profile  on  a  given  body  with  varia- 
tion velocity.   The  corrected  body  profile  radii  decreased 
slightly  with  increase  in  speed.   The  variation  is  considered 
to  be  beyond  the  accuracy  of  the  model  construction  and  the 
full  scale  drawings  so  the  corrected  body  profile  was  plotted 
the  same  for  all  velocities  on  a  given  body.   In  Figures  XVIII 
through  XXII  the  corrected  body  is  plotted  beyond  the  separa- 
tion point  however,  this  in  fact  does  not  exist  since  this 
profile  is  based  upon  modified  potential  flow  theory  which 
does  not  hold  after  separation  occurs. 

The  pressure-velocity  profiles  for  the  bodies  shown 
in  Figures  XIII  through  XVII  are  the  results  of  Computer  Pro- 
gram III.   From  the  sample  of  these  results  given  in  Table  G-II 

U/ 
in  Appendix  G  only  a  slight  variation  in  Cp  and   Uoo  occurs 

with  variation  of  Uoo  or  model  velocity.   This  is  a  result  of 
the  aforementioned  slight  variation  in  body  profile  with  model 
velocity  since  the  output  of  Program  II,  Y  -  NORM,  is  the  in- 
put to  Program  III.   Again  this  slight  variation  is  considered 
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beyond  the  accuracy  of  the  model  drawings  and  certainly  is 
not  detectable  when  plotted  at  the  scales  used  in  Figure  XIII 
through  XVII.  Thus  the  profiles  appear  the  same  for  all  ve- 
locities on  a  given  body. 

Since  neither  depth  of  submergence  nor  hydrostatic 
pressure  enter  into  the  calculations  of  the  corrected  body 
profile  or  pressure-velocity  profiles  none  of  the  figures 
show  a  variation  with  depth. 

In  summary,  the  experimental  results  indicate  that 
separation  appears  to  begin  at  a  critical  angle  of  about 
25  degrees  irrespective  of  body  streamlining.   The  correla- 
tion of  separation  with  speed  and  depth  is  not  clearly  seen 
since  the  speed  and  depth  range  should  be  larger  for  more 
conclusive  results.  The  results  as  pointed  out  in  the  in- 
troduction are  qualitative  in  nature  and  are  not  be  be  con- 
strued as  quantitative  evidence.   Such  evidence  must  be  pro- 
vided in  later  works.   The  primary  recommendation  is  that 
further  studies  be  carried  out  at  greater  speeds  with  im- 
provements in  the  technique  of  observing  the  separation 
point.   St  rat  ford  *s  separation  criterion,  which  is  more  ri- 
gorous than  Buri°s  appears  to  give  good  agreement  with  the 
experimental  observations,  but  it  should  be  further  evaluated 
by  using  experimentally  determined  pressure-velocity  dis- 
tributions. 
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V  -  CONCLUSIONS 

The  conclusions  which  were  arrived  at  in  IV  are  summa- 
rized as  follows: 

1.  The  observed  position  of  separation  was  determined 
within  an  accuracy  of  -  1/2  inch  on  the  body  and  the  body 
angle  at  this  position  was  thereby  determined  to  -  3? 

2.  Initial  separation  occurs  at  body  angles  of  be- 
tween  23°  and  28°  regardless  of  the  amount  of  body  streamlining. 

3.  Buri's  criterion  for  separation  consistently  predicts 
separation  well  forward  of  the  observed  separation  point  and 
thus  does  not  compare  favorably  with  the  observed  results. 

4.  Stratford's  criterion  evaluated  from  the  calcu- 
lated pressure-velocity  distributions  over  the  bodies  gives 
reasonably  close  agreement  with  the  location  of  the  observed 
separation  point. 

5.  Results  of  analyzing  the  separation  position  lo- 
cation show  that  it  moves  aft  with  increase  in  speed  but  this 
may  be  due  to  the  fact  that  at  the  speeds  investigated  the 
Reynolds  numbers  were  in  the  transition  region  between  laminar 
and  turbulent  flow. 

6.  The  theoretical  criteria  evaluated  are  both  in- 
sensitive to  body  speed  changes  and  are  independent  of  the 
depth  of  submergence  of  the  body. 
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7.  The  observed  separation  point  does  not  move  in  a 
consistent  manner  on  the  body  with  a  variation  in  depth  of 
submergence  of  the  body. 

8.  The  corrections  applied  to  the  body  surface  to 
account  for  boundary  layer  effects  showed  little  change  with 
change  in  body  speed  for  a  given  body. 

9.  The  pressure-velocity  profiles  over  the  bodies 
calculated  from  modified  potential  flow  showed  little  change 
with  change  in  body  speed  for  a  given  body. 
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VI  -  RECOMMENDATIONS 

The  recommendations  made  in  IV  are  summarized  as 
follows: 

1.  Further  experimentation  should  employ  a  larger 
model  if  possible  and  an  improved  means  of  observing  the 
flow  about  the  separation  point . 

2.  A  greater  range  of  speeds  and  depths  should  be 
investigated  in  future  work  to  increase  the  reliability  of 
the  data  obtained.   Higher  speeds  consistent  with  the  a- 
bility  to  observe  the  dye-flow  are  necessary  to  insure  that 
flow  is  definitely  in  the  turbulent  region. 

3.  If  higher  speeds  can  not  be  attained  improve- 
ment should  be  made  in  the  method  of  inducing  turbulence 
through  stimulation. 

4.  If  dye  studies  are  continued  a  minimum  of  dye 
should  be  used  to  prevent  dye  outside  the  boundary  layer 
from  obscuring  the  observations. 

5.  Pressure  measurements  should  be  made  in  and 
ahead  of  the  separation  zone  to  experimentally  determine 
the  pressure-velocity  distribution  over  the  bodies  and  also 
serve  as  another  method  of  locating  the  separation  point. 
This  could  be  compared  with  the  results  of  this  work  and 
used  to  re-evaluate  the  separation  criteria  on  the  basis  of 
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the  actual  pressure-velocity  profiles  rather  than  on  calcu- 
lated profiles.   This  appears  to  be  the  next  logical  step 
as  the  initial  determination  of  the  separation  zone  has  been 
made  in  this  work. 

6.   Recommendations  relating  to  drag  force  measure- 
ments are  included  in  Appendix  £. 
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APPENDIX  A 
OFFSETS  AND  OTHER  CHARACTERISTICS  FOR  MODEL  AND  STERN  SECTIONS 


TABLE  A -I  -  OFFSETS 


Body:  Nose  Section   Stern  1  Stern  2 

Station  Spacing:    0.50"       0.50"   0.50" 
Station  0  at:       Nose        Tail    Tail 

Station:       Radius (Inches)  Radius  Radius 


Stern  3  Stern  4  Stern  5 
0.50"    0.50"    0.50" 
Tail    Tail    Tail 

Radius  Radius  Radius 


0 

0.5 

1 

1.5 

2 

3 

4 

5 

6 

7 

7.26 

8 

9 

10 

11 

11.30 

12 

13 

14 

14.78 

15 

16 

17 

18 

19 

20 

21 

22 

22.78 


0.00 

0.43 

0.68 

Uj 

0.83 

0.94 

o 

1.12 

1.27 

o 

1.37 

o 

1.47 

OJ 

II 

1.54 

to 

— — 

z> 

1.62 

H 

1.67 

< 

1.72 

QC 

1.76 

li. 

— 

0 

1.80 

1.83 

111 

1.87 

Ul 

I 

Q_ 

1.90 

in 

1.93 

E 

1.95 

Ul 

1.96 

3: 

1.97 

1.99 

2.00 

0.00 

0,00 

0.00 

0.00 

0.80 

0.65 

0.63 

0.32 

1.05 

0.87 

0.85 

0.47 

1.23 

1.03 

0.98 

0.57 

1.37 

1.16 

1.09 

0.69 

1.58 

1.36 

1.27 

0.87 

1.72 

1.52 

1.41 

1.00 

1.84 

1.66 

1.53 

1.12 

1.93 

1.77 

1.63 

1.23 

1.98 

1.84 

1.72 

1.33 

2.00 

— 

__ 

— 

1.91 

1.78 

1.42 

1.95 

1.84 

1.50 

1.98 

1.89 

1.57 

1.99 

1.93 

1.63 

2.00 

— 

— - 

1.96 

1.69 

1.98 

1.74 

1.99 

1.79 

2.00 

1.83 
1.87 
1.90 
1.92 
1.94 
1.96 
1.98 
1.99 
2.00 
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TABLE  A- I I-  OTHER  CHARACTERISTICS 


Stern 

Stern 

Stern 

Stern 

Body 

Body 

Body 

Body 

Body 


Nose/Tail  Length   Wetted  Surface  Wetted  Surface 
Radius    (inches)   Depths  1  and  2    Depth  3 
(inches) (sq.ft.  ) (sq.ft.  ) 


Nose  0.43 
Cylindrical  Section 

Stern  1  2.00 

2  1.47 

3  0.98 

4  0.99 

5  0.33 


10.5 

19.0 

2.00 

3.63 

5.65 

7.39 

11,39 

31.50 

33.13 

35.15 

36.89 

40.89 


0.726 
1.653 
0.175 
0.245 
0.387 
0.513 
0.741 
2.554 
2.624 
2.766 
2.892 
3.120 


0.640 
1.311 
0,158 
0.219 
0.340 
0.454 
0.664 
2  109 
2,170 
2.291 
2.405 
2.615 


Body  means  total  of  nose,  cylindrical  section,  and  stern  shape. 

Wetted  surfaces  submerged  (depths  1  and  2)  were  calculated 
using  ten-ordinate  Simpson's  rule  with  half-stations  on  body  perimeters 

Spacings  of  1  inch  were  used  on  nose  and  stern  5,  3/4  inch  on 

stern  4,  1/2  inch  on  sterns  2  and  3.   Stern  1  was  calculated 

directly  since  it  was  a  hemisphere. 

Surfaced  wetted  surfaces  (Depth  3)  were  calculated  by  computing 
non-wetted  surface  Using  trapezoidal  rule  on  non-wetted  arc  lengths 
and  subtracting  from  submerged  wetted  surface  values. 

These  calculations  are  not  shown  in  this  report. 
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TABLE      A- III 


GAUSSIAN   STATION, RADIUS 

,    AND  BODY   ANGLE 

Gauss  i: 

an  X 

i      L 

.005 

Body   I 

Body 

2 

Body 

3 

Statioi 

X 

Radius   Angle      X     Radius 

Angle 

s      X     Radius 

Angle 

-.990 

0.17 

0.37 

49.5      0.18 

0.37 

49.4 

0.19 

0.40 

47.3 

-.944 

.028 

0.87 

0.88 

26.0      0.92 

0.91 

25.0 

0.97 

0.94 

23.9 

-.866 

.067 

2.12 

1.30 

14.2      2.23 

1.31 

12.7 

2.36 

1.37 

11.7 

-.756 

.122 

3.85 

1.60 

7.5      4.05 

1.63 

7.1 

4.30 

1.67 

6.3 

-.618 

.191 

6.02 

1.81 

4.5     6.33 

1.84 

4.4 

6.72 

1.87 

3.6 

-.458 

.271 

8.54 

1.96 

2.6      8.98 

1.98 

1.8 

9.52 

1.97 

1.5 

-.282 

.359 

11.31 

2.00 

0.0   11.90 

2.00 

0.0 

12s62 

2.00 

0.0 

-.096 

.452 

14.25 

2.00 

0.0   14.99 

2.00 

0.0 

15.90 

2.00 

0.0 

+  .096 

.548 

17.25 

2.00 

0.0   18.14 

2.00 

0.0 

19.24 

2.00 

0.0 

+  .282 

.641 

20.19 

2.00 

0.0   21.23 

2.00 

0.0 

22.52 

2.00 

0.0 

+  .458 

.729 

22.96 

2.00 

0.0  24.15 

2.00 

0.0 

25,62 

2.00 

0.0 

+  .618 

.809 

25.48 

2.00 

0.0  26.80 

2.00 

0.0 

28.43 

2.00 

0.0 

+  .756 

.878 

27.65 

2.00 

0.0  29.08 

2.00 

0.0 

30.85 

1.94 

5.1 

+  .866 

.933 

29.38 

2.00 

0.0  30.90 

1.78 

13.1 

32.7$ 

1.62 

15.5 

+  .944 

.972 

30.63 

1.63 

35.0  32.21 

1.33 

27.8 

34.18 

1.15 

24.5 

+  .990 

.995 

31.33 

0.78 

67.2  32.95 

0.69 

59.6 

34.96 

0.57 

51.5 

Gauss  ii 

an  X 
i      L 

X        1 

Body 

4 

Body   5 

Statioi 

lad i us 

Angle  X       Radius 

Ang] 

.e 

-.990 

.005 

0.20 

0.41 

47.1      0.21 

0.42 

46.5 

-.944 

.028 

1.02 

0.97 

22.7      1.13 

1.00 

21.6 

A/ore 

-.866 

.007 

2.48 

1.38 

11.5      2.75 

1.42 

10.0 

-.756 

.122 

4.51 

1.67 

5.7      5.00 

1.72 

5.1 

X  is   inches   from 

-.618 

t191 

7.05 

1.88 

3.4      7,81 

1.91 

2.7 

nose   of 

body. 

-.458 

.271 

10.00 

1.99 

1.0  11.08 

2.00 

0.0 

Radius   is 

in 

-.282 

.359 

13y25 

2.00 

0.0  14.69 

2.00 

0.0 

inches. 

-.096 

.452 

16.69 

2.00 

0.0   18.50 

2.00 

0.0 

Angle   is   in 

+  .096 

.548 

20.20 

2.00 

0.0   22.39 

2.00 

0.0 

degrees. 

+  .282 

.641 

23.64 

2.00 

0.0  26.20 

2,00 

0.0 

+  .458 

.729 

26.89 

2.00 

0.0  29.81 

1.98 

1.0 

+  .618 

.809 

29.84 

1.99 

1.0  33.07 

1.85 

4.1 

+  .756 

.878 

32.38 

1.84 

6.2  35.89 

1.57 

8.3 

+  .866 

.933 

34.41 

1.52 

12.0  38.14 

1.18 

12.9 

+  .944 

.972 

35.87 

1.10 

22.0  39.76 

0.73 

20.4 

+  .990 

.995 

36.69 

0.55 

53.0  40.67 

0.28 

36.4 
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APPENDIX  B  -  DESCRIPTION  OF  TOWING  TANK  AND  INSTRUMENTATION 

An  early  description  of  the  Ship  Model  Towing  Tank, 
prior  to  the  installation  of  the  towing  carriage,  is  given  in 
a  paper  entitled  "The  Ship  Model  Towing  Tank  at  M.I.T.,"  Society 
of  Naval  Architects  and  Marine  Engineers,  1953 .   This  paper 
does,  however,  give  an  up  to  date  description  of  the  tank  di- 
mensions and  the  speed  measuring  and  counting  equipment   Since 
the  towing  carriage  was  installed  late  in  1961,  a  brief  de- 
scription will  be  included  here. 

The  carriage  is  constructed  of  aluminum  and  rides  on 
two  hollow  cylindrical  stainless  steel  rails,  one  overhead  and 
one  mounted  on  the  tank's  side  wall,,   The  overhead  rail  is 
simply  for  support  and  guidance  of  the  carriage  which  rides  it 
on  six  fibreglass  wheels,  two  above  and  one  below  the  rail  at 
either  end  of  the  carriage  frame.   The  carriage  is  also  sup- 
ported by  the  side  rail  with  three  wheels  riding  on  the  after 
end  and  a  soft  rubber  drive  wheel  connected  via  a  coupling  to 
an  electric  motor  on  the  forward  end.  The  electric  drive  motor 
is  controlled  by  the  console  at  one  end  of  the  tank  via  cables 
attached  to  pulleys  which  slide  along  a  smaller  rail  mounted 
on  the  tank  wall  above  the  side  rail.   The  control  cables  are 
free  riding  and  consequently  as  the  carriage  moves  from  one 
end  of  the  tank  to  the  other  the  speed  decreases  as  more  cable 
length  is  either  drawn  out  from  one  end  or  piled  up,   A  dif- 
ferential pulley  system  was  installed  which  helped  in  alleviating 
this  situation  to  some  degree. 
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The  carriage  speed  can  be  controlled  either  manually  by 
setting  a  dial  on  the  console  or  automatically  at  a  set  speed 
or  to  position  the  carriage  over  a  self  propelled  model.   The 
automatic  control  was  inoperative  during  the  period  of  the  ex- 
periments included  in  this  thesis,  and  speed  was  controlled 
manually.   The  speed  measuring  apparatus  includes  a  type 
500  Hewlett  Packard  frequency  meter  and  a  type  521  Hewlett 
Packard  counter  described  in  greater  detail  in  the  foremen- 
tioned  reference* 

Figure  B-I  shows  a  schematic  drawing  of  the  model  dye 
system  which  is  self-explanatory  and  figure  B-II  shows  the 
details  of  the  drag  measuring  towing  dynamometer. 
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APPENDIX  C  -  COMPUTER  PROGRAMS 

The  IBM  7090  electronic  computer  located  at  the  M.I.T, 
Computation  Center  was  utilized  as  an  aid  in  reducing  drag 
data  and  in  computing  the  theoretical  portions  of  the  thesis, 

Three  programs  were  written: 
Program  I  -  Reduction  of  drag  data. 

Program  II-  Calculation  of  turbulent  boundary  layer  charac- 
teristics on  the  bodies  tested. 
Program  III-  Calculation  of  the  pressure  velocity  distribu- 
tion on  the  bodies  tested. 

These  programs  were  written  in  7090  FORTRAN.   On  the 
next  several  pages  the  following  are  given  for  each  of  the 
three  programs: 

1.  FORTRAN  program  listing. 

2.  A  list  of  program  symbols  and  their  meanings. 

3.  A  logical  sequence  flow  chart  of  the  program. 
Both  FORTRAN  and  compiled  binary  machine  language  card 
decks  and  listings  are  available  for  these  programs  in  the 
Department  of  Naval  Architecture  and  Marine  Engineering 
computer  program  files. 
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PROGRAM  I  -  LISTING  OF  SYMBOLS 


SYMBOL 

MEANING 

UNITS 

RO 

Mass  density  of  tank  water 

slug/ft 

VIS 

Kinematic  viscosity  of  tank 

water   ft.  /sec 

SLEN(N) 

Body  length 

ft. 

WSUR(N,M) 

Wetted  surface  of  body  at 
particular  depth 

ft.2 

DT(L,M,N) 

Measured  drag 

lbs. 

V' (L,M,N) 

Model  velocity 

kts. 

N 

Body  index 

M 

Depth  index 

L 

Speed  index 

SL 

Lenth  of  strut  submerged 

inches 

DSTR 

Strut  drag 

lb. 

DTC 

Corrected  body  drag 

lb. 

SLR 

Speed  -  length  ratio (V/fh) 

■fW- 

RM 

Reynolds  number  on  SLEN  (N) 

A 

1/2  5  S  V2 

lb. 

CT 

Total  drag  coefficient 

CFSG 

Friftion  drag  coefficient  usj 

•ng 

Shultz-Grunow   formulation 


CFS 


Friction  drag  coefficient  using 
Schoenherr  formulation 


CR 


Residual  drag  coefficient 
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LOGICAL  SEQUENCE  FOR  PROGRAM  1 


START 

I 


READ  f , (RO) 
and  7/,  (VIS) 


READ  LENGTH  OF 
BODY,(SLEN(N)) 


READ  MEASURED 
DRAG,DT(L,M,N) 


READ  VELOCITY 
V(L,M,N) 


SET  BODY  INDEX 

N-l 

=3 


PRINT  BODY 
NUMBER  -  N 


SET  DEPTH  INDEX 
M-l 


®-7 


*r 


PRINT  DEPTH-M 
and  OUTPUT  HEADING 


if  m-2  y^<^> 


STRUT  LENGTH, 
SL-0.0  l 


v 


STRUT  LENGTH, 
SL-6.0  


STRUT  LENGTH, 
SL-14.0 


SET  VELOCITY  INDEX 
L=l 


i 


COMPUTE  STRUT  DRAG, 
(DSTR),  Equation  (27) 


DSTR-0.0 


^COMPUTE  CORR. 
DRAG,  (DTC) 


L 


COMPUTE  V/l'?(SLR) 
and  REYNOLDS  NUMBER, (RN) 


1 


COMPUTE  Ct , (CT) 


[CFS 


COMPUTE  Cf  from  Schultz- 
Grunow  Formula,  (CFSG) 
j Equation 


CFSG-CFS 


COMPUTE  Cf  from  Schoenherr 
Formula  using  CFSG 
Equation 


(IF (ABS.VAL. (CFSG-CFS)  -.OOl) 


COMPUTE  Cr, (CR) 


SPRINT  V(L,M,N)SSLR,DT(L9M,N),DTC, 
CT,CF,CR,  and  RN 


(IS  L-6  >»<^$>»[L-LTl 


(IS  M"3>-<^>HM"M+1 


© 


N-N+ 1  h—<^><IS~N-^>K^E^>  END 
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PROGRAM  I  FOR  CALCULATION  OF  CR  ON  BODY 

FORTRAN 

DIMENSION  DT(6,3,5),V(6,3,5),SLEN(5),WSUR(5,3) 

READ  24, RO,  VIS 

24  FORMAT  (F5.3,  E10.4) 
READ  25, (SLEN(N),N-1,5) 

25  FORMAT  (5F4.2) 

READ  26,  ((WSUR(N,M),  N-l , 5) ,M-1 ,3) 

26  FORMAT  (14F5.3) 

READ  27,  (((DT(L,M,N),  L-1,6)  M-l ,3) ,N-1 , 5) 

27  FORMAT  (18F4.2) 

READ  28,  (((V(L,M,N),  L-1,6),  M=l ,3) ,N=1 ,5) 

28  FORMAT  (18F4.2) 
DD  50  N-l, 5 
PRINT  32,N 

32  FORMAT  (14H1  BODY  NUMBER-T1) 
DP  50  M-l, 3 

PRINT  33, M 

33  FORMAT  (8H0  DEPTH-I1) 
PRINT  34 

34  FORMAT  (111HJ0   VEL.  (KT)      REYNOLDS  NO.     SPEED/LENGTH   TOT.DRA 
1  G    TOT. DRAG (CORR)  CT  CF  CR/) 

IF  (M-2)  10,11,12 

10  SL-14.0 
GO  TO  14 

11  SL-6.0 
GO  TO  14 

12  SL-0.0 

14  DO50  L-1,6 

IF  (M-2)  61,61,60 
64   DSTR-0.0 
GO  To  62 

61  DSTR-(RO*((l.689*V"(LpM,N))**2.0)*((0.064*SL)+1.17))/288.0 

62  DTC-DT(L,M,N)-  DSTR 
SLR-V(L,M,N)/SQRTF(SLEN(N)) 
RN-1.689*V(LJIMJN)*SLEN(N)/VIS 
A-0.5*RO*WSUR(N,M)*(1.689*V(L,M,N))**2.0 
CT-DTC/A 

Xl-0 . 434294*LOGF (RN ) 
CFSG«0.0725/(Xl-2.0)**2.0 

13  X2-0.434294*LOGF(RN*CFSG) 
CFS-0.0585/x2**2.0 
IF(ABSF(CFSG-(0.001E-03))7,7,15 

15  CFSG-CFS 
GO  TO  13 

7  CR-CT-CFS 

PRINT  35,  V(L,M,N),RN,SLR,DT(L,M,N),DTC,CT,CFS,CR 

35  FORMAT  (3H   F4.2,7H        1PE10.3,9H        0PF5.2,1CH 
10PF5.3,8H        OPFS.3,10H  1PE10.3,4H    1PE10.3,5H 
21PE10.3) 

50  CONTINUE 
CALL  EXIT 
END 
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PROGRAM  II  -  LISTING  OF  SYMBOLS 


Symbol 
P 

VIS 

SLEN(N) 

V(L,M,N) 

X(I,N) 

R(I,N) 

AL(I,N) 

N 

M 

L 

I 

RN 

PBLT 

PDT 

BLTO 

BLT 

A 

DLNP 

DLTNP 

DT 


Meaning 

Reciprocal  of  power  in 
velocity  profile 

Kinematic  viscosity  of  tank  water 

Body  length 

Model  velocity 

Distance  from  body  nose 

Radius  of  body  at  X(I,N) 

Cosine  of  angle  between  body  sur- 
face and  body  axis  at  X(I,N) 

Body  index 

Depth  index 

Speed  index 

Location  on  body  index 

Reynolds  number  on  X 

Boundary  layer  thickness  on 
flat  plate  (  6pl) 

Displacement  thickness  on  flat 
plate  (£*  pi) 

Estimate  of  boundary  layer 
thickness  on  body 

Boundary  layer  thickness  on 
body  (6) 

Ratio  of  boundary  layer  thickness 
to  model  radius 

An  pi 


Units 


AZn 


ft.  /sec. 

ft. 

kts. 

inches 

inches 


inches 
inches 

inches 
inches 


Pi 


Displacement  thickness  on  body 
<£*) 


inches 
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PROGRAM  II  (CONT'D) 


Symbol 
OHNP 

OMTNP 

BMT 
RC 


Meaning 


Units 


Pi 

-ft2n 


Pi 


Momentum  thickness  on  body  (©) 

Body  radius  corrected  to 
account  for  displacement 
thickness 

Normalized  radius 


inches 
inches 
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LOGICAL  SEQUENCE  FOR  PROGRAM  II 


START 

m 


© 


READ  l/n,(P);2<(VIS);and 
BODY  LENGTH, (SLEN(N)) 


COMPUTE  ©  ON  FLAT  PLATE 
(BMT),  Equation  (3a) 


READ  VELOCITY 
V(L,M,N) 


I 


SET  S   ON  BODY*  £  ON  FLAT  PLATE 
(BLTO-PBLT) 


SAD  LOCATION  ON  BODY 
X(I,N) 


X 


READ  BODY  RADIUS 
R(I,N) 


COMPUTE  &   ON  BODY 
(BLT ) , Equat  i  on  ( 1 ) 


BLTO-BLT 


^ 


IF  (ABS.VAL.  (BLTO-BLT) -.001) 


READ  COSINE  OF  ANGLE  BETWEEN 
BODY  SURFACE  AND  BODY  AXIS 
AL(N,I) 


I 


SET  BODY  INDEX 

N=l 


D 


or  0s 


COMPUTE  £  *  ON  BODY , 
(DT)„  Equation  (2) 


N-N+l < 


SET  DEPTH  INDEX 
Uzl 


<D 


COMPUTE  0  ON  BODY, 
(BMT),  Equation  (3) 


M-M+l 


-® 


SET  VELOCITY  INDEX 
L-l 


I 


COMPUTE  CORRECTED  BODY  RADIUS, 
(RC),  Equation  (4) 


L-L+l 


fPRINT  OUTPUT  HEADINGS 


SET  LOCATION  INDEX 
1^1 


I-I+l- 


k. 


COMPUTE   REYNOLD'S 
NUMBER,     (RN) 


COMPUTE  NORMALIZED  BODY   RADIUS 
(Y) ,    Equation    (4a) 


<&> 


PRINT  X(I,N),R(I,N),RC,BLT, 
DT,BMT,RN!,and   Y 


(IS   L-3> <yEJfr>^(lS    I-16>^Q>-^-(?) 


COMPUTE  £  ON   FLAT  PLATE, 
(PBLT), Equat ion    (Id  ) 


COMPUTE  <£*   ON   FLAT  PLATE, 
(PDT),    Equation    (2b) 


<JK> — <— (is  m-3; 

<£e^> 
®        <IS  N-5    )^>-^-0 


ES>-^- 


END 
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PROGRAM  II,  CALCULATION  OF  BODY  DISPLACEMENT  THICKNESS , ETC. 

FORTRAN 

DIMENSiaJV(3,3,5),X(l6,5),R(l6,5),SLEN(5),AL(l6,5) 
READ  21,P,VIS,(SLEN(N),N-1,5) 

21  FORMAT  (F3.1,E10.4,5F4.2) 

READ  22, (((V(L,M,M),L-1,3),M-1,3),N-1,5) 

22  FORMAT  (18F4.2) 

READ  23,(<X(I,M),I-1,16),N-1,5) 

23  FORMAT (14F5. 2) 

READ  24,((R(I,N),I-1,16),N«1,5) 

24  FORMAT  (16F4.2) 

READ  25,((AL(I,N),I-1,16),K-1,5) 

25  FORMAT  (14F5.3) 
DO  50  N~l,5 

DO  50  M-1,3 
DO  50  L-1,3 
PRINT  100, N 

100  FORMAT  (14H2  BODY  NUMBER-I1) 
PRINT  110, M 

110  FORMAT  (13H0  BODY  DEPTH- II) 
PRINT  101,V(L,M,N) 

101  FORMAT  (16H0  VELOCITY (KTS)«F4«  2$ 
PRINT  102 

102  FORMAT  (95H0  X(INCH)    RADIUS    CORR.  RADIUS    DELTA    DELTA 
1STAR    THETA     REYNOLDS   NO,   Y-NORM/) 

DO  50  1-1,16 

RN=1.689*V(L,M,N)*X<I,N)/(VIS*12.0) 

PBLT-0.38*X(I,N)/RN**0.2 

PDT-PBLT/(1.0+P) 

BLTO-PBLT 
12  BLT=PBLT/(l.0+(BLTO/(3.0*R(I,N))))**0.8 

IF  (ABSF(BLTO-BLT)-.001)  7,7,10 
10  BLTO-BLT 

GO  TO  12 
7  A-BLT/R(I,N) 

DLNP-PDT/PBLT 

DLTNP-DLNP* (P+l . 0)/ (2. 0*P+1 . 0) 

DT-PDT* (1 . 0+A-SQRTF(l .0+2. 0*A* (1 .0~DLNP)+ (A**2. 0)* (1 . 0-DLTNP) ) )/(A 
1  *DLNP*(l.0+A/3.0)**0o8) 

OMNP-p/ ( (P+l . 0)* (P+2. 0) ) 

OMTNP-OMNP* (P+2. 0)/2 . 0*P+1 . 0) 

BMT-BLT* (OMNP+A*OMTNP/2. 0) 

RC«R(I5N)+AL(I,N)*DT 

Y«RC/(SLEN(N)*6.0) 

PRINT  103,X(I,N),R(I,N),RC,BLT,DT,BMT,RN,Y 

103  FORMAT  (3H    F5.2,6H     F4.2,7H         F5.3,9H       F5.4,7H 
1     F5.4,6H      F5.4,5H     1PE10.3,4H      0PE12.5) 

50  CONTINUE 
CALL  EXIT 
END 
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PROGRAM  III  -  LIST  OF  SYMBOLS 


Symbol 

(D+(J) 
X(I) 


y(i) 

R(l) 

RO 

PATD(M) 

V(L,M,N) 

N 

M 

L 

AL(I) 

A(I) 
B(l) 
G(J,I) 

C(J,I) 
D(J,J) 
SUM(I) 
I(K,I) 
H(I) 

K 


Meaning 


Units 


Gaussian  abscissae,  number  (1-16) 
Normalized  Gaussian  abscissae  based  on  model 
length  of  two,  from  -1  to-*-  1 

Normalized  ordinate  of  body  meridian  profile 

Gaussian  weighting  functions 

Water  density  slug/  ft** 

Static  pressure  (gage)  PsiG 

Absolute  model  velocity  kts. 

Body  number 

Depth 

Speed 

Cosine  of  angle  between  tangent  to  body  profile 

and  x~axis 

K{I)      ratio  of  length  to  diameter 

kjd),  longitudinal  virtual  mass  coefficient 

g(x,t),  equation  of  prolate  spheroid  whose  ends 
coincide  with  ends  of  the  given  body,  and  which 
intersects  body  at  x  -t. 

K  (x,t),  matrix  of  functions  of  body 

K/(x,t)  , Matrix  of  functions  of  prolate  spheroid 

<£(K-K») 

E(n,t), error  functions 

k,  +  £ 
"l  +  k(  (x) 

Iteration  number 
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PROGRAM  III  -  (CONT*D) 


Symbol      Meaning  Units 

SUM(I)      2<En  .  -  Eni) 

nj     *•*■ 

0(1)        ^Eni 

LT(D         V^U*  ,  ratio  of  flow  velocity  over  body  to  free 
stream  velocity. 

P(I)        p/q   pressure  distribution 

VEL(I)       Absolute  flow  velocity  over  body  kts. 

PCC(I)      Gage  pressure  on  body  lb/ft 
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LOGICAL  SEQUENCE  FOR  PROGRAM  III 


Start 


I 


Read  X(I),  R(I),  RO, 
PATD(M),  V(L,M,N) 


Read  N,M,L,  AL(17) 


ris  j  - 16  ") 

»— <no>    <^e^> 


J  -  J+l 


Read  Y(I),AL(I) 


Compute  <£(K-K°) 


P  -  4 


S~\      Compute  E     ,      /"^\ 
(J  /"^  E(1,I)  eq*1  (19)-\3J 


(2  J-^P-P+1 


I  -1 


8 


Compute  H(I)  - 
B(17)  -f  E(l,I) 


1  +2>(\7') 


-H   3 


'3  Wis  1-16   ) 


no> 


I-I+l 


9 
9A 


K  -  1 


J  -  1 


9B 


Compute  :£K(E  •-  Eni) 
feK(J,I)x(E(K,i)-E(K  ,1))) 


x—v  Compute  A.(x)/(A(t) 
[4  h  equation  (17) 


kd 


(is  J  -  16  ) 


^~^\     t 


Compute  K^Cx)  ,  (BCD) 
(  5  V»   equation  (18) 
V_y    Compute  B(17) 


J  -  J-l 


J 


X   Is  I  -  16  ) 


Compute  g(x,t), 
(G(J,I))  eq. (16) 


(Ts  K  -5   ) 

IT 


"n 


Compute  En4.1;L 
(E(K+1,I)) 
eq. (20) 


Compute  K(x,t) , 
(C(J9I)  eq,(l4) 


P  -P  +1 


I  41 


Compute  K' (a,t) , 
(D(J,I))   eq.(l5) 


•® 


K-K+l 
1-1 


(continued  on  next  page) 
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<SH 


Print  N,M,L  and 
output  heading 


€> 


€> 


Compute  Eni 
E(K,I),K-  1..5 


Compute  U(I) 

-(1+1^  +lEni)  cos  <* 

-<1+B(17)+ZE(K,I))AL(I) 


Compute   pj/q 
P(I)-1-U(I)2 


Compute   Vgbs 
VEL(I)-U(I)xV(L,MpN) 


Compute  P  _..,-„ 

PCC(I)-P(I)xR0xV 
2 

+  patdcm: 


M.N) 


2 


Print  X(I),Y(I),U(I),   /Os 

p(i)ppcc(i),vel(i)    ^Viy 
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PROGRAM  III  CALCULATION  OF  PRESSURE  DISTRIBUTION 

FORTRAN 

DIMENSION  X(16),R(  16)  ,Y(16)  ,  AL(16  ,  5)  ,  A(17)  ,B:(17)  ,C(16, 16)  ,D(16, 16) 
1,SUM(16),E(5,16),H(16),0(16),U(16),S(16),P(16),G(16,16),PATD(3), 
2V(3,3,5),VEL(16),PCC(16) 

READ  1,(X(T),R(D,T-1,16) 

1  FORMAT  (6E12.5) 
READ2,RO, (PATD(M) ,M-1,3) 

2  FORMAT (F5 . 3, 3E1 1.4) 
READ3,(((V(L,M,N),L-13),M-1,3),N-1,5) 

3  FORMAT  (18F4.2) 

READ7 , ( (AL(I ,N) ,1-1 , 16) ,N-1 , 5) 
7  FORMAT (14F5. 3) 

4  READ5,NSM9L!)A(17) 

5  FORMAT (311, El 2. 5) 
RE AD6,(Y(l), 1-1,16) 

6  FORMAT (6E12. 5) 
DO  10  1-1,16 
A(r)-SQRTF((X(I)+l.)*(l.-X(I))/(tf(I)**2.)) 

10  CONTINUE 

DO  20  1-1,17 

BCD-(A(l)*LOGF(A(T)+SQRTF(A(I)**  2.-1.  ))-SQRTF(A(I)**2.-l .))/ 

l((A(I)**2.)*SQRTF(A(l)**2.-l.)-A(r)*LOGF(A(r)+SQRTF(A(r)**2,-l.))) 
20  CONTINUE 

DO  30  1-1,16 

DO  30  J-1,16 

G(J!)I)-(Y(T)**2.)*(X(J)+1.)*(1.-X<J))/((X(I)+1.)*<1.-X(T))) 

C(J8T)-R(J)*(Y(J)**2.)/(<(X(J)-X(I))**2.+Y(J)**2.)**1.5) 

D(J,r)-R<J)*G<JfT)/<<(X<J)-X(I))**2.+G<J/I))**1.5> 
30  CONTINUE 

DO  40  1-1,16 

SUM(I)-0.0 

DO  .39  J-1,16 

S(J)-C(J,I)-D(J,I) 

SUM(I)-SUM(I)+S(J) 

39  CONTINUE 
K-l 
E(K9I)-l.-(((t+B(l7))/2.)*SUM(l))-(l.  +  B(l7))/(l.+B(I)) 

H(l)-CE(KJr)+B(l7))/(l.HhB(17)) 

40  CONTINUE 

DO  50  K-l, 4 
DO  50  1-1,16 
SUM(l)-0.0 
DO  49  J-1,16 

8(J)-C(J,r)*(E(K,J)-E(K,l)) 
SUM(r)-  SUM  (I)  *  S(J) 
49  CONTINUE 
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PROGRAM  III  (CONT'D) 


K(Kfl)l)-H(X)*B(K,l)-.5*SUM(r) 
50  CONTINUE 
PRINT  52, N 

52  FORMAT (14H2  BODY  NUMBER-I1) 
PRINT53,M 

53  FORMAT  (8H0  DEPTH-I1) 
PRINT54,V(L,M,N) 

54  FORMAT  (8H0  SPEED- F4. 2) 
PRINT55 

55  FORMAT  (83H0        X(l)         Y(l)         U(T)        P(D 
1         PCC(I)  VEL(X)/) 

DO  60  1-1,16 

0(l)^E(l,I)+E(2,I)+E(3,I)HhE(4,I)+E(5,I) 
U(l)-(l.+B(l7)+0(l))*AL(l,N) 
P(l)-1.-(U(I))**2. 
VEL(I)-U(I)*V(L,M,N)   ' 

PCC(T)-P(I)*(.5*R0*((V(L,M,N)*1.689)**2.))+PATD(M)*144. 
57  PRINT  59,  (X(l)  ,Y(I)  ,U(I)  ,P(r)  ,PCC(I)  ,  VEL(T)) 

59  FORMAT (6E14. 5) 

60  CONTINUE 

GO  TO  (4, 4, 4, 4, 61), N 

61  GO  TO  (4,4,62),M 

62  GO  TO  (4,4,63),L 

63  CALL  EXIT 
END 
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APPENDIX  D  -  EXPERIMENT 
TABLE  D-I  -  EXPERIMENTAL  DRAG  DATA 


Run  Stern  De 

13      1 

1A   ' 

' 

2 

2A 

3 

'■ 

3A 

JL 

4 

2 

4A 

J 

5 

5A 

6 

'' 

6A 

2 

7 

3 

7A 

ii 

8 

8A 

9 

'               .. 

9A   : 

J      3 

10    ] 

I      3 

10A 

11A 

11B 

12 

1 

12A 

3 

13 

2 

13A 

1 

14    ] 

L 

14A 

i 

15 

i 
i 

15A 

2 

16 

m  m 

16A 

, 

17 

17A 

18 

' 

18A   ] 

L       1 

Dial  Speed  Voltmeter  Voltmeter  Voltmeter  Measured 
Depth  Speed  (KTS)   Zero(V)   AVG  (V)   Dlff . (V)  Drag(lbs) 


20 

0.80 

-43.7 

-26.4 

17.3 

0.21 

20 

0.77 

-45.0 

-27.9 

17.1 

0.21 

24 

1.12 

-46.5 

-20.4 

26.1 

0.33 

24 

1.12 

-45.7 

-18.5 

27.2 

0.34 

28 

1.48 

-45.0 

-14.75 

30.25 

0.38 

28 

1.48 

-43.5 

-11.7 

31.8 

0.40 

20 

0.75 

-45.1 

-29.2 

15.9 

0.20 

20 

0.74 

-47.7 

-30.8 

16.9 

0.21 

24 

1.09 

-45/5 

-22.6 

22.9 

0.24 

24 

1.09 

-44.7 

-18.5 

26.2 

0.33 

28 

1.50 

-43.0 

-18.4 

24.6 

0.31 

28 

1.51 

-43.5 

-19.0 

24.5 

0.31 

20 

0.77 

-46.0 

-30.7 

15.3 

0.19 

20 

0.79 

-46.1 

-30.9 

15.2 

0.19 

24 

1.16 

-45.5 

-23.6 

21.9 

0.27 

24 

1.17 

-45.0 

-25.1 

19.9 

0.25 

28 

1.62 

-46.0 

-19.6 

26.4 

0.33 

28 

1.58 

-45.0 

-21.0 

24.0 

0.30 

20 

0.77 

-45.0 

-23.2 

21.8 

0.28 

20 

0.77 

-46.0 

-21.2 

25.0 

0.31 

24 

1.16 

-46.0 

-24.0 

22.0 

0.28 

24 

1.21 

-46.5 

-25.8 

20.7 

0.26 

28 

1.58 

-46.0 

-19.7 

26.3 

0.33 

28 

1.58 

-45.5 

-20.0 

25.5 

0.32 

20 

0.75 

-35.0 

-22.6 

12.4 

0.16 

20 

0.75 

-36.0 

-24.2 

11,8 

0.15 

24 

1.19 

-35.0 

-13.2 

21.8 

0.28 

24 

1.11 

-37.5 

-15.2 

22.3 

0,28 

28 

1.61 

-42.5 

-12.5 

30.0 

0.38 

28 

1.61 

-44.0 

-12.6 

31.4 

0.39 

20 

0.74 

-25.0 

-L5.4 

09  o  6 

0.12 

20 

0.74 

-26.7 

-16.1 

10,6 

0,13 

24 

1.13 

-39,0 

-16d 

22,8 

0,28 

24 

1.16 

-40.5 

-15.2 

25,3 

0.31 

28 

1.56 

-44,0 

-09.2 

34,8 

0.44 

28 

1.53 

=46.0 

-12.6 

33,4 

0.42 
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TABLE  D-I  (CON'T)  EXPERIMENTAL  DATA 


Dial   Speed  Voltmeter  Voltmeter  Voltmeter  Measured 
Run  Stern  Depth   Speed  (KTS)  Zero(V)    AVG  (V)   Diff . (V)   Drag(lbs) 

0.19 
0.20 
0.32 
0.36 
0.39 
0.32 
0.16 
0.14 
0.21 
0.21 
0.23 
0.23 
0.28 
0.26 
0.26 
0.26 
0.27 
0.26 
0.25 
0.29 
0.31 
0.31 
28     1.65   -40.0      -11.2      28-8       0.35 

0.35 
0.20 
0.21 
0.26 
0.25 


19   2     1 

19A 

/i 

20 

20A 

21 

" 

21A 

1 

22 

2 

22A  2 

23 

23A 

24 

' 

24A 

2 

25 

3 

25A 

* 

26 

26A 

27 

r 

27A  2     3 

28  4     3 

28A 

1      ' 

29 

29A 

30 

>r 

30A 

3 

31 

2 

31A 

■  i' 

32  4 

32A 

33 

" 

33A 

2 

34 

1 

34A 

» 

35 

35A 

36 

r 

36A  4     1 

37   5     1 

37A 

J 

38 

38A 

20 

0.75 

-44.0 

20 

0.73 

-45.0 

24 

1.14 

-46.0 

24 

1.16 

-46.5 

28 

1.61 

-44.5 

28 

1.56 

-43.5 

20 

0.80 

-24.0 

20 

0.73 

-24.5 

24 

1.15 

-28.0 

24 

1.14 

-26.0 

28 

1.67 

-25.0 

28 

1.60 

-24.5 

20 

0.73 

-37.0 

20 

0.75 

-33 . 0 

24 

1.18 

-32.0 

24 

1 .  18 

-30.0 

28 

1.53 

-29.0 

28 

1.50 

-26.0 

20 

0.70 

-40.0 

20 

0.70 

-42.5 

24 

1.00 

-45.0 

24 

1.25 

-42.5 

28 

1.65 

-40.0 

28 

1.67 

-40.0 

20 

0.78 

-35.0 

20 

0.78 

-34.5 

24 

1.22 

-31.5 

24 

1.23 

-32.0 

28 

1.63 

-32.0 

28 

1.58 

-30.0 

20 

0.77 

-43.0 

20 

0.76 

-48.0 

24 

1.20 

-47.0 

24 

1.18 

-40.5 

28 

1.58 

-41.0 

28 

1.58 

-44.0 

20 

0.78 

-37.0 

20 

0.80 

-44 . 0 

24 

1.21 

-30.5 

24 

1.18 

-34.5 

28.6 

15.4 

29.4 

15.6 

20.0 

26.0 

17.9 

28.6 

13.2 

31.3 

17.9 

25.6 

09.3 

14.7 

11.3 

13.2 

09.0 

19.0 

07.5 

18.5 

04.6 

20.4 

04.0 

20.5 

13.2 

23.8 

10.3 

22.7 

09.4 

22.6 

07.6 

22.4 

05.5 

23.5 

03.8 

22.8 

08.6 

21.4 

17.5 

25.0 

19.2 

25.8 

16.2 

26.3 

11.2 

28.8 

11.2 

28.8 

17.2 

17.8 

15.7 

18.8 

08.8 

22.7 

10.0 

22.0 

07.4 

24.6 

05.2 

24.8 

19.7 

23.3 

24.1 

23.9 

15.2 

31.8 

12.3 

28.2 

05.4 

35.6 

07.7 

36.3 

18.1 

18.9 

20.9 

23.1 

05.0 

25.5 

06.2 

28.3 

0, 
0< 
0, 
0, 
0, 
0, 
0, 
0, 


29 
30 
26 
28 
39 
34 
45 
45 


0.21 
0.27 
0.30 
0.34 
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TABLE  D- 1 (CONT'D)  EXPERIMENTAL  DATA 


Dial 

Speed 

Voltmeter 

Voltmeter 

Voltmeter  Measured 

Run 

Stern  Depth 

Speed 

(KTS) 

Zero  (V) 

AVG 

(V) 

Diff . 

(V) 

Drag (lbs) 

39 

J 

t 

28 

1.67 

-37.5 

-01.1 

■ 

36.4 

0.46 

39A 

i 

28 

1.80 

-43.0 

-07,3 

35.7 

0.44 

40 

2 

20 

0.80 

-45.0 

-20.1 

24.9 

0.29 

40A 

5 

I 

20 

0.80 

-46.0 

-20.3 

25.7 

0.30 

41 

24 

1.28 

-46.5 

-15.3 

31.2 

0.38 

41A 

24 

1.28 

-46.0 

-13.5 

32.5 

0.40 

42 

■ 

r 

28 

1.70 

-46.0 

-  8.8 

37.2 

0,47 

42A 

4 

I 

28 

1.65 

-46 . 0 

-  8.4 

37.6 

0.47 

43 

r 

J 

20 

0.75 

-43.5 

-18.1 

25.4 

0.30 

43A 

1 

k 

20 

0.75 

-39.5 

-17.3 

22.2 

0.26 

44 

24 

1.25 

-39.0 

-14.4 

24.6 

0.29 

44A 

24 

1.28 

-38.5 

-13.1 

25.4 

0.30 

45 

< 

1 

t 

28 

1.63 

-37.0 

-07,4 

29,6 

0.35 

45A 

5   : 

I 

28 

1.60 

-36,0 

-07.0 

29.0 

0.35 

Depth  1  -  18"  to  Keel  (4.5*Dia.) 
Depth  2  -  10"  to  Keel  (2.5*Dia,) 
Depth  3  -  3,6**to  Keel  (0.9*Dia.) 

Runs  -  1  -  21A  Conducted  24  March  1962 
22  -  45A  Conducted  25  March  1962 


Water  Temp  -  51°  F 


1.386  x  10"5  ft2/sec 


S    -  1.939  lb  -  sec2/ft 
Dial  speed  is  setting  on  manual  speed  control  knob 


80 


TABLE  D-II  - 

SEPARATION  POINT 

DATA 

Run     Stern        Run 

Separation 

Body 

Classificat 

ion 

Point 

Angle 

1        3         LD 

5.0 

33° 

1A 

1          LD 

5.0 

33 

2 

MD 

5,0 

33 

2A 

MD 

5.0 

33 

3 

HD 

5,0 

33 

3A 

HD 

5.0 

33 

4 

LM 

4.75 

28 

4A 

LM 

4.75 

28 

5 

MM 

4.75 

28 

5A 

MM 

4.75 

28 

6 

HM 

4.75 

28 

6A 

HM 

5.0 

33 

7 

LS 

4.5 

23 

7A 

LS 

4,5 

23 

8 

MS 

4.5 

23 

8A 

MS 

4.5 

23 

9        < 

r            HS 

4.75 

28 

9A      : 

3         HS 

4.5 

23 

10 

L         L3 

0.8 

24 

10       I 

LS 

0  =  8 

24 

11A 

MS 

1.25 

42 

11B 

MS 

1.25 

42 

12 

HS 

1.3 

42 

12A 

HS 

1.35 

42 

13 

LM 

0.75 

24 

13A 

LM 

0.75 

24 

14 

MM 

1.0 

32 

14A 

MM 

1.1 

32 

15 

HM 

1.1 

32 

15A 

HM 

1 

.1-1.2 

32 

16 

LD 

0 

.8-09 

24 

16A 

LD 

0.8 

24 

17 

MD 

0.8 

24 

17A 

MD 

0 

,8-0,9 

24 

18       i 

'          HD 

1.3 

42 

18A 

L          HD 

1.3 

42 

19 

2                      LD 

2.75 

30 

19A      i 

LD 

2.5 

25 

20 

MD 

3.0 

36 

20A 

MD 

3.1 

36 

21 

HD 

? 

- 

21A 

HD 

3.0 

36 
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TABLE  D-II    (cont'd) 


Run     Stern          Run 

Separation 

Body 

Classification 

Point 

Angle 

22 

L1C 

3.0" 

36° 

22A 

LM 

2.6 

25 

23 

MM 

2.6 

30 

23A 

MM 

2.8 

30 

24 

HM 

? 

«, 

24A 

HM 

2.75 

30 

25 

LS 

2.9 

33 

25A 

LS 

2.8 

30 

26 

MS 

3.1 

36 

26A 

MS 

3.2 

42 

27       1 

HS 

O  •  <3 

42 

27A      2            HS 

3.25 

42 

28       4            LS 

7.0 

38 

28A 

I             LS 

7.0 

38 

29 

MS 

7.1 

44 

29A 

MS 

7.1 

44 

30 

HS 

7.3 

49 

30A 

HS 

7.3 

49 

31 

LM 

7.1 

44 

31A 

LM 

7.1 

44 

32 

MM 

7.2 

49 

32A 

MM 

7.2 

49 

33 

HM 

7.2 

49 

33A 

HM 

7.2 

49 

34 

LD 

6.8 

28 

34A 

LD 

6.7 

28 

35 

MD 

6.8 

28 

35A 

MD 

6.9 

.30 

36 

»             HD 

7.0 

38 

36A 

1             HD 

? 

_ 

37       ! 

5             LD 

10.5? 

23 

37A 

LD 

? 

=> 

38 

MD 

10.5 

23 

38A 

MD 

10.5 

23 

39 

HD 

10.5 

$ 

39A 

HD 

10.5 

23 

40 

LM 

11.0 

30 

40A 

LM 

11.0 

30 

41 

MM 

? 

_ 

41A 

MM 

? 

= 

42 

HM 

10.9 

30 

42A 

HM 

? 

=_ 

43 

LS 

? 

- 

43A 

LS 

? 

- 

44 

MS 

11.5 

60 

44A 

MS 

11.5 

60 

45 

HS 

11.5 

60 

45A      J 

3             HS 

11.5 

60 
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TABLE  D-II  (Cont'd) 
NOTE: 

Run  classification  code  -  1st  letter  designates  speed 

L  -  low 
Id  -  medium 
H  -  high 

2nd  letter  designates  depth 

D  -  depth  1  (deep) 
M  -  depth  2  (medium) 
S  =  depth  3  (surfaced) 

Separation  point  is  measured  in  inches  from  end  of  cylindrical 
section.   Marks  on  photographs  are  \   inch  apart. 

Body  angle  is  angle  between  body  axis  and  body  surface  at 
separation  point. 
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APPENDIX  E  -  DISCUSSION  OF  DRAG  FORCE  MEASUREMENTS 

The  experimental  data  from  Table  D-I  was  reduced  to  co- 
efficient form  using  the  following  procedure. 
(Drag  measured)  -  (Strut  drag)  -  (Total  corrected  drag) 
where  strut  drag  was  calculated  by  the  method  that  follows. 
Measurement  of  a  strut  tare  drag  was  attempted  but  due  to 
either  the  end  effects  on  the  strut  because  of  the  lack  of 
a  dummy  body  or  faulty  instrumentation  these  values  were 
of  a  larger  magnitude  than  the  values  of  measured  drag  of 
the  body  plus  strut  and  thus  were  not  valid.   An  engineer- 
ing  value  of  strut  drag  composed  of  profile  drag,  spray 
drag,  wave  drag,  and  interference  drag  with  the  body  was 
calculated  as  follows: 
Profile  Drag: 

t/c  -    .8  -   .246   for  strut  used. 
3755 

Reynold  Number  on  Strut  *  3  x  104 

from  figure  126,  page  247  of  (9) 

Cd  -  0.08 

Drag  -  q  Cd  (.8  x  SL^  -  q  (.064)  SL 
144  144 

where  SL  -  strut  length  submerged  (inches) 
2  /,  _\2 


q-lP  Vk   (1.689) 


Spray  Drag:  _j 

Froude  number  on  strut  **"  4  x  10 

t  -   .8  -  .246 
2x    3.25 

from  figure   25,    page   10-13  of    \l6] 
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Cdx  -  .24  /t\2=.058  x  -  1.625 

Drag  -  q  Cdx(X2)  -  q   (.1535) 

144 

Wave  Drag: 

from  figure  24,  page  10-13  of  (lOj 

Cdt^  1.4 

Drag  -  q  Cdt  t2  -_q_  (.895) 

144 

Interaction  Drag: 

from  figure  23,  page  8-10  of  [lOj 

Cdt  =   .75  x  t/c  =   .19 

Drag  -  Cd  q  t2  -  q   (.121) 
at       tu 


Total  Strut  Drag: 


Strut  Drag  -  1  j>   VK    (1.689)2 


.064SL  +1,17 


(27) 


2  144     L  J 

The  total  corrected  drag  was  converted  to  coefficient 
form  by  the  following  relationship: 

C+  -    Drag  (28) 

t   2 

1/2  ?   VK   (1.689)2  (g) 
where  S  -  wetted  surface  of  the  model. 
The  coefficient  of  frictional  resistance,  C^ ,  was  nejt  calcu- 
lated using  Schoenherr's  relationship  that: 


Cf=  0.0585 


log1Q  (RftCf)J 


-2 


-  0.0585  o  (29) 

r.4343   loge  (Rn  Cf)J 

In  order  to  use  the  electronic  computer  to  evaluate  (29) 


the  Schultz-Grunow  formulation 
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Cp  -  0.0725  )log10  Rn  -2,o7 
-  Q.0725 


-2 
2  <30) 


[.4343  logeRn  -2.0 J 
which  yields  Cf  values  greater  than  that  of  (29)  in  the  range 
of  interest  in  this  problem  was  calculated  and  used  as  a  first 
value  to  solve  (29)  which  is  implicit  in  Cf. 

The  residual  resistance  coefficient,  Cr,  was  then 
calculated  from: 

Cr  -  Ct  -  Cf  (31) 

It  was  hoped  that  perhaps  some  correlation  between  the 

change  of  location  of  separation  on  the  body  and  change  in  Cr 

for  various  runs  could  be  found  if  our  drag  measurements  were 

good. 

An  IBM  7090  comuter  program,  Program  I  was  written  to 

solve  equations  (27)  through  (31).  The  details  of  this  program 
are  given  in  Appendix  C.   Sample  outputs  of  this  program  are 
included  in  Table  E-I  at  the  end  of  this  appendix.   These  values 
which  are  typical  of  the  values  obtained  on  all  the  bodies  are  much 
too  large  for  the  underwater  shapes  tested.   At  submerged  depth 
1  the  values  of  Cr  and  C^  are  over  ten  times  as  large  as  would 
be  expected  by  comparison  to  Series  58  data. 

At  first  strut  drag  was  suspected  but  in  the  surfaced  con- 
dition where  strut  drag  does  not  enter  into  the  calculation  the 
residual  drag  values  are  much  too  large  also.   The  bow  waves 
generated  in  the  surfaced  condition  are  shown  in  Figure  IX.   It 
is  not  possible  for  the  wave  making  resistance  created  in  this 
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condition  to  account  for  such  large  values  of  C  .   Besides 
being  too  large  in  magnitude  the  values  are  not  consistent 
from  body  to  body  with  the  highly  streamlined  body  some- 
times showing  a  higher  C  value  than  a  less  streamlined  body 
for  similar  speed  and  depth  of  submergence.   Plots  of  C 
versus  both  Reynolds  number  and  speed-length  ratio  were 
prepared  but  these  were  completely  inconsistent  and  unsatis- 
factory when  compared  to  similar  plots  of  similar  bodies  in 
the  literature. 

At  the  time  of  this  writing  the  cause  of  the  excess 
drag  is  unknown.   It  is  speculated  that  the  dynamic  response 
of  the  drag  measuring  system  is  causing  a  false  reading  of 
drag  due  to  some  phenomena  that  does  not  appear  during  sta- 
tic calibration.   The  calibration  curves  used  in  this  work 
are  shown  in  Figure  £-1 .   It  is  also  of  interest  that  a^ 
nother  group  using  another  submerged  model  experienced  the 
same  difficulty,  i.e.,  too  much  drag. 

If  it  is  desired  to  make  further  drag  measurements 
on  the  model  it  is  recommended  that  the  drag  force  dyna- 
mometer be  mounted  at  the  model  rather  than  in  the  strut. 
This  will  eliminate  the  problem  relating  to  strut  drag  with 
the  exception  of  the  interference  drag  between  the  strut  and 
body.   It  is  also  recommended  that  the  dynamic  response  of 
the  entire  drag  measuring  system  used  during  the  tests  be 
given  a  thorough  investigation.   The  authors  realize  that 
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the  new  towing  tank  and  dynamometer  equipment  was  put  to 
use  on  an  emergency  basis,  the  construction  of  such  appara- 
tus being  completed,  and  before  the  carriage  and  instrumen- 
tation were  properly  calibrated,  accepted,  and  put  into 
regular  use.   Upon  completion  of  all  of  these  steps  it 
should  be  feasible  to  carry  out  successful  drag  measure- 
ments on  submerged  bodies  at  the  new  M.I.T.  Towing  Tank. 
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Figure  EtI-  Dynamometer  Calibration 
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TABLE  E-I  -  SAMPLES  OF  COMPUTER  OUTPUT,    PROGRAM  I 


BCCY    MJI"fcER=l 


U6PTH: 

=  1 

VEL.  1 

;kt 

0.74 

C.74 

1.13 

1.16 

1.56 

1.53 

DEPTH: 

-1 

VEL.I 

;kt 

0.75 

0.75 

l.iy 

l.  ll 

1.61 

1.61 

DEPTH= 

=  3 

VEL.I 

;kt 

0.77 

0.77 

1.16 

1.21 

1.58 

1.58 

REYNOLDS  NO, 
2.3  72E  o5 
2.372L  05 
3.6226  05 
3.718E  0  5 
5. COOt  05 
4.904E  0  5 


REYNOLDS    NO 

2. 404E  05 

2.404E  05 

3.8  UE  0  5 

3.558t  05 

5.160E  05 

5.160t  05 


SPI 

EEO/LENGTH 

TOT. DRAG 

TOT.DRAu(CORR) 

C.46 

0.120 

C.098 

C.46 

0.  130 

C.  108 

C.70 

0.280 

0.229 

0.72 

0.310 

0.257 

C.96 

0.440 

0.343 

C.94 

0.420 

0.327 

SPEED/Lt 

NGTH 

TOT. DRAG 

C.46 

C.  160 

0.46 

C.  150 

0.73 

0.280 

0.68 

0.280 

C.99 

0.38C 

0.99 

0.390 

TOT.ORAO(CORR) 
0.143 
0.  133 
0.238 
0.243 
0.303 
0.313 


DEPTH 


REYNOLDS  MO. 

SPEED/LENGTH 

TOT. DRAG 

TOT.DRA 

2.468t  05 

0.47 

0, 

.280 

0.280 

2.468E  05 

0.47 

0, 

.310 

0.310 

3.718E  05 

C.72 

c, 

.280 

0.280 

3.878E  05 

0.75 

0. 

.260 

0.260 

5.C64E  05 

C.97 

o. 

.330 

0.330 

5.C64E  05 

0.97 

0, 

.320 

0.320 

CT 

CF 

CR 

2.541E- 

-02 

5.909E- 

-03 

1.950E- 

-02 

2.799E- 

-02 

5.909E- 

■03 

2.208E- 

-02 

=  1    2.543E- 

-02 

5.400E- 

-03 

2.003E- 

-02 

2.700E- 

-02 

5.371E- 

-03 

2.163E- 

-02 

1.998E- 

-02 

5.053E- 

-03 

1.493E- 

-02 

1.978E- 

-02 

5.073E- 

-03 

1.471E- 

-02 

DEPTH  =  2 


CT 
3.604E-02 
3.353E-02 
2.377E-02 
2.795E-02 
1.653E-02 
1.707E-02 


CF 
5.892E-03 
5.892E-03 
5.342E-03 
5.421E-C3 
5.021E-03 
5.021E-03 


CR 
3.015E-02 
2.763E-02 
1.842E-02 
2.253E-02 
1.151E-02 
1.205E-02 


DEPTH  =  3 


CT 
7. 797E-02 
8.632E-02 
3.436E-02 
2.932E-02 
2.182E-02 
2.  116E-02 


CF 
5.858E-03 
5.858E-03 
5.371E-03 
5.324E-03 
5.040E-03 
5.040E-03 


CR 
7.211E-02 
8.047E-02 
2.898E-02 
2.400E-02 
1.678E-02 
1.612E-02 
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BODY  NUKRER=3 


U6PTb=l 


TABLE  E-I   (CON'T) 


VEL.  (KT, 

0.8C 

0.77 

1.12 

1.12 

l.4b 

1.48 


RLYNLLDS    NO, 

2.8  56E  0  5 
2.749E  05 

3.9  99E  0  5 
3.9  99E  0  5 
5.284b  05 
5.284t  05 


SPEED/LENGTH 

TOT. DRAG 

TOT.DRAG(CGRR) 

0.47 

0.210 

0.  185 

C.45 

0.210 

0.  186 

C.65 

0.330 

0.280 

C.65 

0.340 

0.290 

C.86 

0.38C 

0.293 

C.86 

0.400 

0.313 

l)EPTh  =  2 


VEL.UT, 

0.75 

0.74 

1.09 

1.C9 

1.5C 

1.51 


REYNOLDS  NO, 

2.678E  05 

2.642E  05 

3.892t  05 

3.892E  05 

5.356t  05 

5.392E  05 


PEED/LENGTH 

TOT. DRAG 

C.44 

0.200 

0.43 

0.210 

0.64 

0.290 

C.64 

0.330 

0.88 

0.310 

C.88 

0.310 

TOT. DRAG (CORR) 
0.183 
0.194 
0.255 
0.295 
0.243 
0.242 


DEPTHS 


0.77 
0.79 
1.16 
1.17 
1.62 
1.58 


KT. 

REYNOLDS    NO. 

SPEED/LENGTH 

TOT. DRAG 

TOT.DRA 

2.749E    05 

C.45 

0 

.  190 

0.190 

2.82U    05 

C.46 

0 

.  190 

0.190 

4-.-1426    05 

C.68 

0 

.270 

0.270 

4.178t     05 

C.68 

0 

.250 

0.250 

5.7&4E    05 

0.95 

0 

.330 

0.330 

5. 6416    05 

C.92 

0 

.300 

0.300 

CT 

CF 

Ck 

3.7716- 

■02 

5.6786- 

-03 

3.2036- 

-02 

4.1116- 

-02 

5.7246- 

-03 

3.5396- 

-02 

-J3EPTH 

=  1            2.9206- 

■02 

5.2906- 

-03 

2.3916- 

-02 

3.0246- 

-02 

5.2906- 

-03 

2.4956- 

-02 

1.7496- 

-02 

4.9976- 

-03 

1.2496- 

-02 

1.8686- 

-02 

4.9976- 

-03 

1.3696- 

-02 

CT 

CF 

CR 

4.258E-02 

5.756E-03 

3.6826-02 

4.623E-02 

5.7736-03 

4.0466-02 

DEPTH  ■  2 

2.801E-02 

5.320E-03 

2.2696-02 

3.241E-02 

5.3206-03 

2.7096-02 

1.4116-02 

4.9836-03 

9.1266-03 

1.387E-02 

4.9766-03 

8.8946-03 

DEPTH  =3 


CT 
5.0586-02 
4.6056-02 
3.1676-02 
2.8826-02 
1.9856-02 
1.8976-02 


CF 
5.7246-03 
5.6936-02 
5.2526-03 
5.2436-03 
4.9066-03 
4.9316-03 


CR 
4.4856-02 
4.2356-02 
2.6426-02 
2.3586-02 
1.4946-02 
1.4036-02 
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BODY  NUKBER=5 


DEPTH=1 


TABLE  E-I  (CON'T.) 


VEL.tKT, 

0.78 

0.8C 

1.21 

1.18 

1.67 

1.8C 


REYNOLDS  NO 
3.241E  05 
3.324t  05 
5.028E  05 
4.903E  05 
6.940t  05 
7.4806  05 


PEEO/LENGTH 

TOT. DRAG 

TOT.DRAG(CCRR) 

0.42 

0.210 

0.  186 

0.43 

0.270 

0.245 

0.66 

0.300 

0.242 

0.64 

0.340 

0.285 

C.90 

0.460 

0.349 

C.97 

0.440 

0.311 

0EPTH=2 


VEL.tKT 

0.8C 

0.8C 

1.28 

1.28 

1.70 

1.65 


REYNOLDS  NO 

3.324E 

05 

3.324E 

05 

5.319E 

05 

5.319E 

05 

7.064E 

05 

6.857E 

05 

SPEED/LENGTH 

TOT. DRAG 

0.43 

0.290 

0.43 

0.300 

C.69 

0.380 

C.69 

0.400 

0.92 

0.470 

C.89 

0.470 

TOT.DRAG(CORR) 
0.271 
0.281 
0.331 
0.351 
0.384 
0.389 


DEPTh=3 


VEL.tKT. 

0.75 

0.75 

1.25 

1.28 

1.63 

1.60 


REYNOLDS  NO. 
3.1l7t  05 
3.1171:  05 
5.194L  05 
5.319E  05 
6.773E  05 
6.649E  05 


SPEED/LENGTH 
0.41 
0.41 
C.68 
0.69 
C.88 
0.87 


TOT. DRAG 
0.300 
0.260 
0.290 
C.300 
0.350 
0.350 


TOT.DRAG(CORR) 
0.300 
0.260 
0.290 
0.3C0 
0.350 
0.350 


DEPTH 


CT 
3.540E-02 
4.429E-02 
1.915E-02 
2.370E-02 
1.452E-02 
1.U4E-02 


CF 
5.528E-03 
5.498E-03 
5.047E-03 
5.073E-03 
4.731E-03 
4.662E-03 


CR 
2.987E-02 
3.879E-02 
1.410E-02 
1.863E-02 
9.785E-03 
6.478E-03 


-DEPTH 


CT 
4.905E-02 
5.086E-02 
2.342E-02 
2.483E-02 
1.539E-02 
1.655E-02 


CF 
5.498E-03 
5.498E-03 
4.990E-03 
4.990E-03 
4.714E-03 
4.742E-03 


CR 
4.355E-02 
4.537E-02 
1.843E-02 
1.984E-02 
1.067E-02 
1.181E-02 


DEPTH 


CT 
7.374E-02 
6.391E-02 
2.566E-02 
2.532E-02 
1.821E-02 
1.890E-02 


CF 
5.574E-03 
5.574E-03 
5.014E-03 
4.990E-03 
4.754E-03 
4.771E-0.' 


CR 
6.817E-02 
5.834E-02 
2.065E-02 
2.033E-02 
1.346E-02 
1.413E-02 
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APPENDIX  F  -  SAMPLE  CALCULATIONS 

Calculation  of  Hydrostatic  Pressure. 

The  hydrostatic  pressure  to  the  keel  of  the  bodies  at 
the  three  depths  of  submergence  is  desired  for  use  in  com- 
puter Program  III  to  compute  the  gage  pressure  on  the  body 
from  the  following  relationship; 


Pgage  -  /£\Q  +  Ps  (144) 


where  P  -  dynamic  pressure  (ps£) 

%   -  1/2  <^U  2   (psf  ) 

P$  ■  hydrostatic  pressure  (psi) 

Ps  -  <^9  h 
1728" 

where  P  -  1.935  lb. sec2   at  S1°F 

9  -  32.174  ft/sec2 

h  -  depth  to  keel  (inches) 
Using  the  foregoing  relation  the  following  results  were 
obtained: 


Depth 

h (inches) 

Ps(psi) 

1 

18 

.6484 

2 

10  -,. 

.3602 

3 

3.6 

.1297 

?3 


Calculation  of  Point  of  Separation  from  Stratford  Criterion 

A  typical  calculation  of  the  location  of  the  separation 
point  according  to  Stratford's  criterion  is  given  below  for  body 
#3  at  velocity  of  0.75  knots. 

From  the  form  of  (21)  separation  can  occur  only  when 
Cp  ,  as  plotted  in  Figure  F-II,  is  positive  since  negative 
values  of  Cp  will  result  in  imaginary  numbers. 

At  X  -  34.8"  X  -  34.8  -0.99 
Z       35.15 

from  Figure  F-I;   U  °  0.75 

Uoo 

from  Table  (b-I;   Reynolds  number  -  2.65  x  105  at  X-34.8" 

thus  Rs  -  .75  (2.65  x  105)  -  1.99  x  105 

from  Figure  F-II;  Cp  -  0.05  and  dCp  -  0.50 

dx 

n  -  log  10  Rs  -  5.30  from  (22) 

o 
since,      dp    ?o;/6a    .73        from   (24) 

dx2 
thus  the  left  hand  side  of  (21)  becomes 


5.30-2 


i 


(2x.05)^~4"    (34.8  x  .50)z  -  .625 

while  the  right  hand  side  of  (21)  becomes 

(1.06) (.73) (1.99  x  10~5xl06)    1  -  .912.   Thus  the 

criterion  is  not  satisfied  at  X  ■■  34.8"  since  the  RHS  is  much 

larger  than  LHS.   Next  try 

X  -  34.9*  X  =  .995 
L 
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from  Figure  F-II;    Cp  -  O.10  and  dCp  -  0.96 


^ 


thus  the  LHS  of  (21)  becomes 

(2  x.10)  4  (34.9x.96)  -  1.46  which  now  makes  the 
RHS  smaller  than  LHS  so  the  point  of  separation  lies  between 
X  -  34.8"  and  X  -  34.9. " 

Calculations  similar  to  the  above  were  used  to  cal- 
culate the  separation  points  on  the  other  bodies  at  the  va- 
rious speeds.   During  these  calculations  it  was  noted  that  the 
criteria  was  insensitive  to  speed  changes  and  that  the  same 
point  of  separation  occurred  for  all  values  of  speed.   The 
points  of  separation  according  to  this  criterion  are  marked 
by  small  circles  on  Figure  F-II. 

Calculation  of  Point  of  Separation  from  Buri's  Criterion 

A  typical  calculation  of  the  location  of  the  separation 
point  according  to  Buri's  criterion  is  given  below  for  body  #3 
at  a  velocity  of  1.12  knots* 

From  the  form  of  equation  (26)  separation  will  occur 

only  when  dV U\  as  determined  from  Figure  F-I  is  negative. 

VUoo/ 

<l 

at  X  -  .906; 

Z 

d/  U\  -  -.115  and  U  -  1.103;  from  Figure  F-I 
^Uoo '  Uoo 

and  Q   -  .091;   from  Table  G-I 
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thus  the  left  hand  side  of  (26)  becomes 

-(.115)(1.103)(.  906)(091)  -  -.0105  >    —.006 

Thus  the  point  of  separation  lies  slightly  forward  of  X  -  .906 

L 
but  in  the  region  of  where 

d/_U_\  is  still  negative.   Separation  point  is  at 

{UooJ 


d(r) 

X  -  .900  or  X  -  31.7  inches  which  is  as  close  as 
Z 
the  point  can  be  calculated  within  the  accuracy  of  the  plot  in 

Figure  F-I. 

Calculations  similar  to  the  above  were  used  to  calculate 

the  separation  points  on  other  bodies  at  various  speeds.   During 

these  calculations  it  was  noticed  that  the  criterion  was  insensitive 

to  speed  changes  and  that  the  same  point  of  separation  occurred 

for  all  values  of  speed.   The  points  of  separation  according  to 

this  criterion  are  marked  by  small  circles  on  Figure  F-I. 

NOTE: 

Figures  F-I  and  F-II  on  the  next  two  pages  were  plotted 
from  values  obtained  from  the  outputs  of  Programs  II  and  III 
as  aids  to  evaluate  (21)  and  (26). 
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FIGURE  F-I   -  PLOT  OF  U     va      X 
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FIGURE  F-II   -  Plot   «f  Cj>  ys   X   (distance   from  body  nose' 
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APPENDIX  G   -    SAMPLES   OF  COMPUTER    OUTPUT.    PROGRAMS gand    HI 


TABLE  Q-I    -   SAMPLES   OF  COMPUTER   OUTPUT,    PROGRAM  H 


BODY  NUMHEK=3 


8O0Y  f)EPTH=l 


VELOCITYIKTS)=l. I? 


X( INCH) 

RADUIS 

CORK.  RADIUS 

DELTA 

DELTA-STAR 

THETA 

REYNOLDS  NO. 

Y-NORM/ 

0.19 

0.40 

0.401 

.0154 

.0019 

.0015 

2.161E  03 

0.22825E-01 

0.97 

0.94 

0.946 

.0564 

.0068 

.0056 

1.103E  04 

0.53822E-01 

2.36 

1.37 

1.383 

.1141 

.0135 

.0114 

2.684E  04" 

"O.  78683E-01 

4.30 

1.67 

1.691 

.1832 

.0214 

.0184 

4.891E  04 

0.96202E-01 

6.72 

1.87 

1.900 

.2599 

.0298 

.0263 

7.643E  04 

0. 10806E-00 

9.5  3 

1.97 

2.003 

.3408 

.0383 

.0348 

1.084E  05 

0.  U424E-00 

12.62 

2.00 

2.047 

.4225 

.0465 

.0437 

1.435E  05 

0. 11641E-00 

15.90 

2.00 

2.054 

.5032 

.0542 

.0526 

1.808E  05 

0.  U685E-00 

19.  24 

2.00 

2.061 

.5806 

.0613 

.06  14 

2.188E  05 

0.  11725E-00 

22.52 

2.00 

2.066 

.6528 

.0677 

.0697 

2.561E  05 

0.  11762E-00 

25.62. 

2.00  • 

2.073 

.7181 

.0734 

.0773 

2.914E  05 

0.11794E-00 

26.43 

2.00 

2.078 

.7751 

.0782 

.0841 

3.234E  05 

0.11821E-00 

30.85 

1.94 

2.C21 

.8206 

.0814 

.0899 

3.509E  05 

0.11497E-00 

32.79 

1.62 

1.697 

•5.42  7 

.0802 

.0947 

3.729E  05 

0.96544E-01 

34.  17 

1.15 

1.217 

.8326^ 

.07  34 

.0985" 

3V886E  05 

0.69213E-01 

34.96 

U.57 

0.605 

.7529 

.0560 

.  1022 

3.976E  05 

0.34408E-01 

HOLSY    NUML(EK  =  3 
eODY    LEPIH=2 
VELCCITY(MS)  =  1.51 


XI  INCH) 

KADUIS 

COKR.  RADIUS 

DtLTA 

DELTA-STAR 

THETA 

REYNOLDS  NO. 

Y-NORM/ 

0.  19 

0.40 

0.401 

.0145 

.0018 

.0014 

2.914E 

03 

0.22821E-01 

0.97 

0.94 

0.946 

.0531 

.0064 

.0053 

1.487E 

04 

0.53R03E-01 

2.36 

1.37 

1.383 

.  1076 

.0128 

.0107 

3.619E 

04 

0.78642E-01 

4.30 

1.67 

1.690 

.172U 

.0202 

.017*3 

6.594E 

04 

0.96138E-01 

6.72 

1.87 

1.898 

.2453 

.0282 

.0248 

1.030E 

05 

0.10797E-00 

9.53 

1.97 

2.006 

.3218 

•  0  3"'.4 

.0328 

1.461E 

05 

0.11413E-00 

12.62 

2.00 

2.044 

.3991 

.0442 

.0411 

1.935E 

05 

0.11628E-00 

15.90 

2.00 

2.052 

.4756 

.0516 

.0495 

2.438E 

C5 

0.  11670E-00 

19.24 

2.00 

2.058 

.5490 

.0585 

.0578 

2.950E 

05 

0.  11709E-00 

22.52 

2.00 

2.065 

.6175 

.0646 

.0656 

3.453E 

05 

0.11744E-00 

25.62 

2.00 

2.070 

.6795 

.0701 

.0728 

3.929E 

05 

0.  11775E-00 

28.43 

2.00 

2.075 

.7333 

.0747 

.0792 

4.360E 

05 

0.  11801E-00 

30.85 

1.94 

2.018 

.7771 

.0779 

.Of  46 

4.731E 

05 

0.U477E-00 

32.79 

1.62 

1.694 

.7987 

.0769 

.0891 

5.028E 

05 

0.96362E-01 

34.  17 

1.15 

1.214 

.7905 

.0706 

.0927 

5.240E 

05 

0.69068E-01 

34^96 

0.57 

0.604 

.7173 

.0542 

.096  1  ' 

5.361E 

05 

0.34343E10f 
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BOCY  NUMBEK=1 
BODY  DEPTh=2 
VEL0CITY(KTS)=0.75 


X(  INCH) 

RAOUIS 

CORR.  RADIUS 

DELTA 

DELTA-STAR 

THETa 

REYNOLDS  NO. 

Y-NORM/ 

0.17 

0.37 

0.  37  I 

.0152 

.0019 

.0015 

1.295E 

03" 

0.23524E-01 

0.87 

o.ao 

0.88o 

.0559 

.0067 

.0055 

6.626E 

03 

0.56149E-01 

2.12 

1.30 
1.60 

1.309 
1.621 

.1134 
.1815 

.0134 
.0211 

.0113 
.0182 

1.615E 

04 

0.82974E-01 

.3.85 

2.932E 

04 

0.10272E-00 

6.02 

1.81 

1.839 

.2576 

.0295 

.0261 

4.585E 

04 

0.  11656E-00 

8.53 

1.96 

1.998 

.3379 

.0380 

.0346 

6.497E 

04 

0.  12661E-00 

11.31 

2.00 

2.046 

.4195 

.0462 

.0433 

3.614E 

04 

0.  12967E-00 

14.25 

2.00 

2.054 

.4997 

.0539 

.0522 

1.085E 

05 

0.  13016E-00 

17.  24 

2.00 

2.061 

.5765 

.0610 

.0609 

1.313E 

05 

0.13061E-00 

20.19 

2.00 

2.067 

.6484 

.0674' 

.0692 

1.538E 

05 

0.  13101E-00 

22.96 

2.00 

2.073 

.7131 

.0729 

.0767 

1.749E 

05 

0.13137E-00 

25.48 

2.00 

2.078 

.7699 

.0777 

.0835 

1.94  IE 

05 

0.  13167E-00 

27.65 

2.00 

2.082 

.8173 

.0816 

.0892 

2.106E 

05 

0.13191E-00 

29.38 

2.00 

2.C85 

.8543 

.0846 

.0937 

2.238E 

05 

0.13210E-00 

30.63 

1.63 

1.697 

.862(1 

.0818 

.0972 

2.333E 

05 

0.  10754E-00 

31.33 

0.70 

0.805 

.7922 

.0637 

.  1005 

2.386E 

05 

0.50993E-01 

BODY  NU^bEK=3 


BODY  DEPTh=2 


VEL0CITY(KTS)=0.75 


X(  INCH) 

RAOUIS 

CCRR.  RADIUS 

DELTA 

DELTA-STAR 

THETA 

REYNOLDS  NO. 

Y-NORM/ 

0.19 

0.40 

0.401 

.0167 

.0020 

.0016 

1.447E 

03 

0.22831E-01 

0.97 

0.94 

0.94  7 

.0610 

.0073 

.0060 

7.388E 

03 

0.53850E-01 

2.36 

1.37 

1.384 

.1234 

.014fc 

.0123 

1.797E 

04 

0.78741E-01 

4.3U 

1.67 

1.693 

.1980 

.0230 

.0199 

3.275E 

04 

0.96293E-01 

6.72 

1.87 

1.902 

.2808 

.0320 

.0285 

5.118E 

04 

0.10819E-00 

9.53 

1.97 

2.011 

.3679 

.0410 

.0378 

7.258E 

04 

0.  11439E-00 

12.62 

2.00 

2.050 

.4559 

.0497 

.0474 

9.612C 

04 

0.U659E-00 

15.90 

2.00 

2.058 

.5426 

.0579 

.0570 

1.21  IE 

05 

0.11706E-00 

19.24 

2. CO 

2.065 

.6256 

.0654 

~;0565 

T.465E 

05 

Tj.ll748E-0tT 

22.52 

2.00 

2.072 

.7030 

.0721 

.0756 

1.715E 

05 

0.  U787E-00 

25.62 

2.00 

2.G73 

.7730 

.0780 

.0839 

1.951E 

05 

0.  11820E-00 

28.43 

2.00 

2.  Ofi  3 

.8340 

.0830 

.0912 

2. 165E 

05 

0.  U849E-00 

30.85 

1.94 

2.026 

.882o 

.0363 

.0975 

2.350E 

05 

0.  11524E-00 

32.79 

1.62 

1.702 

.9051 

.0848 

.  1027 

2.497E 

05 

0.96797E-01 

34.17 

1.15 

1.220 

.892  3 

.0772 

.  1069 

2.603E 

05 

7)y&9*l4"F^0~V 

34.96 

0.57 

0.606 

.8028 

.0586 

.1110 

2.663E 

05 

0.34498E-01 

BODY  NUMBERS 


BODY  CEPTh=3 


VEL0CITY(KTS)=0.75 


X(  INCH) 

RADUIS 

CORR.  RADIUS 

DELTA 

DELTA-STAR 

THETA 

REYNOLDS  NO. 

Y-NORM/ 

0.22 

0.42 

0.422 

.0187 

.0023 

.0018 

1.676E 

03 

0.20604E-01 

1.13 

1.00 

1.008 

.0639 

.0082 

.0068 

8.606E 

03 

0.49250E-01 

2.75 

1.42 

1.436 

.1392 

.0163 

.0139 

2.094E 

04 

0.70191E-01 

5.00 

1.72 

1.746 

.2228 

Y02-57 

."0225 

~3.'308E 

04 

"0V8T316E-01 

7.81 

1.91 

1.946 

.3155 

.0356 

.0322 

5.948E 

04 

0.95092E-01 

11.08 

2.00 

2.046 

.4130 

.0456 

.0426 

8.439E 

04 

0.99979E-01 

14.69 

2.00 

2.055 

.5113 

.0550 

.0535 

I. 119E 

05 

0.  10044E-00 

18.50 

2.00 

2.064 

.6076 

.0638 

.0645 

1.409E 

05 

0.10087E-00 

22.39 

2.00 

2.072 

.7000 

.0718 

.0752 

1.705E 

05 

0.  10126E-00 

26.20 

2.00 

2.079 

.7857 

.0790 

.0854 

1.995E 

05 

0.10161E-00 

29.81 

1.98 

2.065 

.8626 

.0851 

.0948 

2.270E 

05 

0.  10093E-00 

33.08 

1.85 

1.938 

.9239 

.0887 

.1033 

2.519E 

05 

0.94742F-01 

35.89 

1.57 

1.657 

.9614 

.0882 

.  1106 

2.734E 

05 

0.81000E-01 

38.14 

1.18 

1.260 

.9655 

.0824 

.  1169 

2.905E 

05 

0.61602E-01 

39.75 

0.73 

0.795 

.9154 

.0692 

.  1225 

3.027E 

05 

0.38848E-01 

40.67 

3.28 

0.316 

.7423 

.0443 

.1297 

3.098E 

05 

0.1T470F-01 

100 


tiOfiY    MUMHEK=1 


TABLE  G-II  -  SAMPLES  OF  COMPUTER  OUTPUT,    PROGRAM  III 


CEPTH=2 
SPEtl)  =  0.75 

X{  I  ) 

-0.98940E  UO 

-0.94458E  00 

-0.86563E  00 

-0.75540E  00 

-0.61788E  00 

-0.4b802E-00 

-0.28160E-00 

-0.95013E-01 

0.9S013E-01 

0.28160E-U0 

C.4b8C2E-00 

0.61788E  00 

0.7b540E  00 

0.86563E  00 

0.94458E  00 

0.98940E  00 


Y(I  ) 
0.23524E-C1 
0.b6149E-Cl 

0.82974E-C1 
0.  LC2  72E-CO 
0.  1  16S6E-C0 
0.12661E-C0 
0.  l?y67E-00 
0.  13016F-C0 
0.130olE-C0 
0.13101E-CO 
0.  13137E-C0 
0.131o7E-C0 
0.13191G-00 
0.  13210E-C0 
0.  107i>4E-C0 
0.b0993E-Cl 


U(  I  ) 
0.6G402E  00 
0.98608E  00 
0.74893E  00 
0.10463E  01 
0.103b8E  01 
0.104  13E  01 
0.10222E  01 
0.10106E  01 
0.1009CE  01 
0.10099E  01 
0.10136E  01 
0.10200E  01 
0.105b2E  01 
0.  12431E  01 
0.11424E  01 
0.55872E  00 


P  (  I  > 

0.55907E  00 
0.27646E-01 
0.43910E-00 
■0.94643E-01 
•0.  728UE-01 
-0.84395E-01 
0.44932E-01 
-0.21251E-QI 
•0.  17995E-01 
-0. 19858E-01 
■0.27446E-01 
-0.40472E-01 
■0.  11343E-00 
-0.54533E  00 
■0.30503E-00 
0.68783E    00 


PCC( I > 
0.52739E  02 
0.51912E  02 
0.52552E  02 
0.51722E  02 
0.51756E  02 
0.51738E  02 
0.131799E  02 
0.51836E  02 
0.51341E  02 
0.51838E  02 
0.51826E  02 
0.51806E  02 
0.51692E  02 
0.51020E  02 
0.51394E  02 
0.52939E  02 


VEUl)/ 
0.49802E-00 
0.73956E  00 
0.56170E  00 
0.78469E  00 
0.77682E  00 
0.78101E  00 
0.76666E  00 
0.75793E  00 
0.75672E  00 
0.75741E  00 
0.76022E  00 
0.76503E  00 
0.79140E  00 
0.93233E  00 
0.85678E  00 
0.41904E-00 


HODY  NUMt'ER=3 
LEPTH=2 
SPEED=0. 75 


X(  1  ) 

-0.98940E  00 

-0.944b8fc  00 

-0.86563E  00 

-0.75540E  00 

•0.61788E  00 

-0.45802E-00 

-0.28160E-00 

-0.95013E-01 

0.95'013E-01 

0.28160E-00 

0.45602E-00 

0.61788E  00 

0.75b40E  00 

0.86563E  00 

0.94458E  00 

0.98940E  00 


Y(  I  ) 
0.22831E-C1 
0.b38bUE-Cl 
0.7S741E-C1 
0.9C293E-C1 
0.  1C819E-C0 
C.  1  1439E-C0 
r.  U659E-CC 
0.  117U6E-C0 
0.  11748E-C0 
C.  U787E-C0 
C.  11720E-C0 
0.  I  1849E-C0 
0.  1  lb24E-C0 
0.96797E-C1 
0.69414E-C1 
C.  34498E-C1 


U(  I  ) 

0.70115E  00 

0.10042F  01 

0.10542E  01 

0.10473E  01 

0.10410E  01 

0.1C281E  01 

0.10166F  01 

0.10088E  01 

0.10081E  01 

0.10123E  01 

0.100G2E  01 

0.10435E  01 

0.11069E  01 

0.10737E  01 

0.10298E  01 

0.77644E  CO 


P(  I) 

0.50839E  00 
-C.85C96E-02 
-0.  1 1 1  2i>E-00 
-0.96933E-01 
-0.83723E-01 
-0.56907E-01 
-0. 33406E-01 
-0.  17689E-01 
-0.  16231F-01 
-0.24782E-01 
-0.  16448E-01 
-0. 88884E-01 
-0.22517F-00 
-0.  15292E-00 
-C.60430E-01 

0. 39714E-00 


PCC(  I  ) 

0.52660E  02 
0.51856E  02 
0.51696E  02 
0.51718E  02 
0.51739E  02 
0.51780E  02 
0.51817E  02 
0.51841E  02 
0.bl844E  02 
0.bl830E  02 
0.51643E  02 
0.bl731E  02 
0.51518E  02 
0.bl63lE  02 
0.51775E  02 
0.b2487E  02 


VEL(  I  )/ 

0.52586E  00 

0.75318E  00 

0.79062E  00 

0.7Bbi>lE  00 

0.78077E  00 

0.77104E  00 

0.76242E  00 

0.75660E  00 

0.7b606E  00 

0.7b924E  00 

0.75614E  00 

0.78262E  00 

0.83016E  00 

0.80531E  00 

0.77233E  00 

0.b8233E  00 
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BODY    NUMHEK=j 


OEPTH=l 

SPEED=l. 12 

X(  I  ) 

Y(  1  ) 

U(  I  ) 

P(I) 

PCC(  I  ] 

i 

VEL(  I  )/ 

-0.98940E 

GO 

0. 

22825E- 

-01 

0. 

70131E 

00 

G. 

508166 

00 

0.951336 

02 

0. 

78547E 

00 

-0.94458E 

00 

0. 

53B22E- 

-01 

0. 

10043E 

01 

-0. 

861266- 

-02 

0.933406 

02 

0. 

112486 

01 

-0.86563E 

GO 

0. 

76663E- 

-CI 

0. 

10542E 

01 

-0. 

111266- 

-00 

0.92984E 

02 

0. 

U807E 

01 

-0.75540E 

00 

0. 

96202E- 

-01 

0. 

10473E 

01 

-0. 

96911E- 

-01 

0.93033E 

02 

0. 

11730E 

01 

-0.61788E 

00 

0. 

10806E- 

-co 

0. 

10409E 

01 

-0. 

83437E- 

-01 

0.93080E 

02 

0. 

11658E 

01 

-0.458026- 

-00 

0. 

11424E- 

-co 

0. 

102806 

01 

-0. 

56869E- 

-01 

0.93172E 

02 

0. 

115146 

01 

-0.281606- 

•GO 

0. 

11641E- 

-co 

0. 

1C165E 

01 

-0. 

33273E- 

-01 

0.93254E 

02 

0. 

11385E 

01 

-0.95013E- 

-CI 

0. 

11685E- 

-co 

0. 

10087E 

01 

-0. 

174216- 

-01 

0.93309E 

02 

0. 

11297E 

01 

0.95013E- 

-Gl 

0. 

11725E- 

-co 

0. 

10O79E 

01 

-0. 

157896- 

-01 

0.93315E 

02 

0. 

U288E 

01 

0.281606- 

-00 

0. 

1176  2E- 

-co 

0. 

10096E 

01 

-0. 

19289E- 

-01 

C.93303E 

02 

0. 

11307E 

01 

0.A5802E- 

-GO 

0. 

11  7  94E- 

-cc 

0. 

1C149E 

01 

-G. 

29985E- 

-01 

0.93266E 

02 

0. 

11367E 

01 

0.61788E 

GO 

0. 

11821E- 

-co 

0. 

103996 

01 

-0. 

81449E- 

-01 

0.93087E 

02 

0. 

11647E 

01 

0.755406 

GC 

0. 

11497E- 

-cc 

0. 

11064E 

01 

-0. 

22402E- 

-00 

0.92592E 

02 

0. 

12391E 

01 

0.86563E 

00 

0. 

96544E- 

-01 

0. 

10732E 

01 

-0. 

151696- 

-00 

0.928436 

02 

0. 

12019E 

01 

0.94458E 

00 

0. 

69213E- 

-Gl 

0. 

10289E 

01 

-0. 

587316- 

-01 

0.93166E 

02 

0. 

115246 

01 

0.989406 

00 

0. 

34408E- 

-CI 

0. 

77605E 

00 

0. 

39775E- 

-00 

0.947506 

02 

0. 

86917E 

00 

BODY  NUM6EK=3 

DEPTH=2 

SPEED=1.51 

x(  i ) 

Y(  I  ) 

U(l  ) 

P(  I) 

PCC( I  ) 

V6L(  I  )/ 

-0.98940E 

00 

0. 

.22821E- 

-CI 

0. 

70143E 

00 

0. 

50799E 

00 

0.55072E 

02 

0, 

,  105926 

01 

-0.94458E 

00 

0. 

.538G3E- 

-01 

0. 

.10043E 

01 

-0. 

,871556- 

-02 

0.518146 

02 

0. 

,151666 

01 

-0.86563E 

00 

0. 

.78642E- 

-01 

0. 

,  10542E 

01 

-0. 

,  11127E- 

-00 

0.51167E 

02 

0. 

159186 

01 

-0.75540E 

00 

0. 

,961386- 

-01 

0. 

10473E 

01 

-0. 

•  96905E- 

-01 

0.51258E 

0  2 

0. 

,  158156 

01 

-0.61788E 

00 

0. 

.1C797E- 

-co 

0. 

,10409E 

01 

-0. 

•  83397E- 

-01 

0.51343E 

02 

0. 

,  15717E 

01 

-0.45802E- 

-00 

0. 

,  U413E- 

-co 

0. 

,  1G280E 

01 

-0. 

.56824E- 

-01 

0.51510E 

02 

0. 

,  15523E 

01 

-0.281606- 

-00 

0. 

.116286- 

-co 

0. 

,  10165E 

01 

-0. 

,  33229E- 

-01 

0.516596 

02 

0. 

,  15349E 

01 

-0.950136- 

-0  1 

0. 

.116706- 

-GO 

0. 

,  100866 

01 

-0. 

,  173296- 

-01 

0.517606 

02 

0. 

,152306 

01 

0.95013E- 

-01 

0. 

.11709E- 

-co 

0. 

100796 

01 

-0, 

,  15799E- 

-01 

0.51769E 

02 

0. 

152196 

01 

0.28160E- 

-00 

0. 

.  117446- 

-co 

0. 

i  100956 

01 

-0. 

,  19173E- 

-01 

0.51748E 

02 

0. 

,152446 

01 

0.45802E- 

-00 

0. 

.11775E- 

-co 

0. 

.10148E 

01 

-0. 

.29892E- 

-01 

0.51680E 

02 

0. 

,  15324E 

01 

0.61788E 

00 

0. 

.118G16- 

-co 

0. 

,  103986 

01 

-0. 

,811806- 

-01 

0.51357E 

02 

0. 

,157016 

01 

0.75540E 

00 

0. 

.  11477E- 

-co 

0. 

110626 

Gl 

-0. 

223586- 

-00 

0.50459E 

02 

0. 

16703E 

01 

0.86563E 

00 

0. 

■  96362E- 

-CI 

0. 

,107296 

01 

-b. 

,  151056- 

-00 

0.509166 

02 

0. 

,  16200E 

01 

0.94458E 

00 

0. 

.69068E- 

-01 

0. 

,102846 

01 

-0. 

,  57650E- 

-01 

0.515G5E 

02 

0. 

155296 

01 

0.93940E 

00 

0. 

,  34  34  3E- 

-01 

0. 

,775806 

00 

0. 

,  398  136- 

-00 

0.54379E 

02 

0. 

11715E 

01 

BODY  NUMBEk=5 

CEPTh-=3 

SPEED=0.75 

X(  I  ) 

Y(  1  ) 

um 

P(  I) 

pcc( i : 

1 

VEL(  I  )/ 

-0.98940E 

00 

0. 

206U46- 

-Cl 

0. 

70371E 

00 

0. 

50479E 

00 

0. 

194626 

02 

0. 

,527786 

00 

-0.94458E 

00 

0. 

49250E- 

-CI 

0. 

102856 

01 

-G. 

57780E- 

-01 

0. 

18587E 

02 

0. 

771366 

00 

-0.8&563E 

00 

0. 

70191E- 

-01 

0. 

104696 

01 

-0. 

960356- 

-01 

0. 

185276 

02 

0. 

785196 

00 

-0.755406 

GO 

0. 

85316E- 

-CI 

0. 

104356 

01 

-0. 

889186- 

-01 

0. 

18538E 

02 

0. 

782636 

00 

-0.61788E 

00 

0. 

,950926- 

-CI 

0. 

103516 

01 

-0. 

713886- 

-01 

0. 

18566E 

02 

0. 

,77o3lE 

00 

-0.45802E- 

-00 

0. 

,999796- 

-CI 

0. 

10267E 

01 

-0. 

54101E- 

-01 

0. 

185936 

02 

0. 

77002E 

00 

-0.28160E- 

-00 

0. 

,  100446- 

-CO 

0. 

10093E 

01 

-0. 

186116- 

-01 

0. 

18648E 

02 

0. 

75b95E 

00 

-0.95G13E- 

-01 

0. 

1C087E- 

-CO 

0. 

10071E 

01 

-0. 

142366- 

-01 

0. 

18655E 

02 

0. 

Ib-iili: 

00 

0.95013E- 

-0  1 

0. 

,  10126E- 

-CO 

0. 

10076E 

01 

-0. 

152126- 

-01 

0. 

18653E 

0  2 

0. 

,755681 

00 

0.28160E- 

-00 

0. 

.10161E- 

-co 

0. 

10133E 

01 

-0. 

,  26690E- 

-01 

0. 

18635E 

02 

0. 

75  J94E 

00 

0.45802E- 

-00 

0. 

,  10093E- 

-co 

0. 

10  309E 

01 

-C. 

62655E- 

-01 

0. 

18D796 

02 

0. 

773146 

00 

0.61788E 

00 

0. 

.94742E- 

-CI 

0. 

10432E 

01 

-0. 

.883686- 

-01 

0. 

18539E 

02 

0. 

78244E 

00 

0.75540E 

00 

0. 

.81000E- 

-CI 

0. 

103466 

01 

-0. 

,  70363E- 

-01 

0. 

185676 

02 

0. 

,775946 

00 

0.86563E 

00 

0, 

.61602E- 

-CI 

0. 

,  10135E 

01 

-0. 

,272486- 

-01 

0. 

18634E 

02 

0. 

,760151 

00 

0.944586 

00 

0. 

.38848E- 

-01 

0. 

958106 

00 

0. 

,820386- 

-01 

0. 

188046 

02 

0. 

,718586 

00 

0.98S*406 

00 

c. 

,  1543CE- 

-CI 

0. 

771  18E 

CO 

0. 

,404986- 

-00 

0. 

193076 

02 

0. 

,57851t 

CO 
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